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ABSTRACT 


Constant  load  creep  data  were  obtained  for  four  plastics:   polyethyl- 
ine,  polypropylene,  nylon,  and  polycarbonate.   Viscoelastic  responses  of 
these  materials  proved  to  be  non-linear  in  both  time  and  stress.   An 
empirical  creep  equation,  based  on  a  time  dependent  retardation  rate 
constant,  was  developed  to  characterize  the  creep  behavior  and  good  cor- 
relation was  obtained  for  various  stress  levels.   By  assuming  the  de- 
veloped relationship  to  be  an  equation  of  state,  creep  behavior  of  simple 
composites  of  plastics  was  predicted  with  good  results  when  viscoelastic 
responses  of  the  components  were  similar.   The  importance  of  stress-time 
history  is  demonstrated  by  using  the  creep  equation  to  estimate  mechanical 
behavior  in  other  modes  of  testing.   The  concept  of  a  retardation  spectrum 
and  its  application  to  the  materials  studied  is  considered. 
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NOTATION 

"y  -  shear  stress 

V  -  shear  strain 

7?    -  viscosity  coefficient 

y   -  shear  strain  rate 

Gt    -    shear  modulus 

E.  -  tensile  modulus 

t  -  time 

£  -  engineering  strain 

€.  -  strain  rate 

0~  -  engineering  stress 

C  -  stress  rate 

Cq  -   initial  stress 

01 1  -   stress  at  time  "  Lt  " 
3(t)-  creep  compliance 

O  -   current  time 

-£  -  natural  logarithm  base 

la.  -  glass  transition  temperature 

yt   -  free  energy  of  deformation 

Co-  zero  time  strain  (instantaneous  elastic  strain) 

C^-  delayed  elastic  strain 
IBt-(t)-  relaxation  modulus 

€^.  -  non-recoverable  flow  strain 
b   -  deformation 
£  -  deformation  rate 
PE  -  polyethylene 


PP  -  polypropylene 

PC  -  polycarbonate 

NY  -  nylon 

7"  -  retardation  time  constant 

t  -  rate  constant  for  delayed  strain  response 

<2Vt>)  -  rate  constant  distribution  function 

L(TJ  -  retardation  time  constant  distribution  function 

C.   -  delayed  strain  fraction 

A( -  creep  parameter  representing  zero-time  strain 

A?-  creep  parameter  representing  steady-state  creep  rate 

A3-  creep  parameter  representing  magnitude  of  delayed 
elastic  strain  component 

C,   -  time  dependent  rate  constant  in  empirical  creep 
equation 

Q  -  time  exponent  in  empirical  creep  equation 

K^t^i  -  stress  coefficients  in  generalized  creep  equation 

NjM^R  -  stress  exponents  in  generalized  creep  equation 
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1.    Introduction. 

Composite  materials  with  high  strength-to-weight  ratios  show  great 
promise.  [3]   Composites  may  be  looked  upon  as  intimate  mechanical  combina- 
tions of  two  or  more  materials,  but  differ  from  multiphase  materials  such 
as  metallic  alloys  in  that  they  are  not  homogeneous  on  a  macroscopic  scale. 
Composite  materials  are  not  new.   Concrete  and  filled  resins,  for  example, 
have  been  successfully  employed  for  years,  and  currently  a  bewildering 
number  of  combinations  of  many  thousands  of  materials  is  available.   Of 
necessity,  however,  the  work  with  composites  to  date  has  been  directed 
towards  development  of  materials  to  meet  specific  technological  require- 
ments, and  very  little  research  has  been  carried  out  in  the  area  of 
characterization  of  the  behavior  of  composites.  [1,  2] 

Physical  properties  often  constitute  the  starting  point  for  the 
materials  selection  process,  and  so  characterizing  mechanical  behavior 
of  composites  is  an  important  first  step  in  a  fundamental  research  study. 
Of  particular  interest  is  viscoelastic  behavior,  frequently  manifested  by 
the  component  materials  of  the  composite  structure.   Deviations  from 
true  elasticity  are  negligible  for  engineering  purposes  for  materials 
such  as  metals  and  glasses,  but  in  polymeric  systems,  which  are  extensively 
employed  in  composites,  behavior  is  dominated  by  viscoelastic  phenomena. 
[2]  Many  conventional  composite  materials  have  mechanical  properties 
that  exhibit  appreciable  time  dependence  through  creep  and  relaxation  pro- 
cesses. [4]   The  stress  analysis  becomes  quite  complicated,  but  it  must 
be  performed  if  a  viscoelastic  composite  is  to  be  utilized  in  an  effici- 
ently designed  structure  giving  satisfactory  performance  in  service 
environments.  [5] 

In  this  investigation  the  main  purpose  was  to  characterize  the  creep 
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behavior  of  individual  viscoelastic  materials  by  analytic  equations,  and 
then  obtain  suitable  correlations  so  that  the  behavior  of  composites  of 
these  materials  could  be  predicted. .  In  addition,  these  analytic  expressions 
were  used  to  estimate  performance  of  materials  in  modes  of  testing  other 
than  creep,  such  as  stress  relaxation  and  constant  strain  rate  uniaxial 
tension  tests. 

For  the  experiment,  four  different  thermoplastics  (polyethylene,  poly- 
propylene, nylon,  and  polycarbonate)  were  chosen,  and  composites  were 
formed  as  layered  sandwich  structures.   A  simple  composite  was  selected 
for  this  basic  study  to  eliminate  as  many  variables  as  possible,  since 
stress  analysis  of  polymeric  materials  is  inherently  complex,  even  if 
linear  viscoelasticity  is  assumed.   In  all  experiments  the  specimens  were 
subjected  to  uniaxial  tension,  and  it  has  been  shown  [6]  for  this  type  of 
loading  that  changing  the  type  and  amount  of  component  interfacial  con- 
tact has  negligible  effect  on  composite  behavior.   This  investigation  was 
carried  out  in  a  laboratory  with  controlled  temperature  of  23+1°  C. 
and  relative  humidity  50  +  2°L   in  accordance  with  ASTM  Standards.  [7] 
Laboratory  set-up  and  equipment  is  reported  in  Appendix  G. 
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2.    Linear  Viscoelasticity. 

For  the  ideal  elastic  solid,  stress  is  proportional  to  strain  when 
a  load  is  applied  in  accordance  with  Hooke's  Law.   It  is  also  independent 
of  the  rate  of  strain.   For  example,  at  stresses  below  the  yield  strength, 
the  behavior  of  most  metals  is  given  by  Hooke's  Law: 

(2.i)  7-=  GT 

where   /  =  applied  shear  stress 
T  -   shear  strain 
Q  =  modulus  of  rigidity 

For  the  ideal  viscous  liquid,  stress  is  proportional  to  the  rate  of 

strain,  but  independent  of  the  strain  itself,  as  stated  by  Newton's 

Viscosity  Law.  [8]   At  low  shear  stresses  most  liquids  and  gases  obey 

Newton's  Law: 

(2.2)       J-    =  fff 

where  J=   applied  shear  stress 
0   =  shear  strain  rate 
T\   =  Newtonian  viscosity 
coefficient 

The  laws  and  equations  above  represent  idealizations.   If  sufficiently 
precise  measurements  are  made  it  is  found  that  no  real  solid  is  truly 
Hookean  and  that  no  real  liquid  is  exactly  Newtonian  in  mechanical  re- 
sponse. [9]   To  represent  such  deviations  an  equation  can  be  written  in 
the  form: 

(2.3)  Y    =    BTt^ 

For  the  ideal  elastic  solid,  />t  equals  zero  and  the  equation  becomes: 

(2.4)  Y    =    57 
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where  the  constant  D  is  the  reciprocal  of  the  shear  modulus,  /Q  .   For 
a  perfectly  viscous  liquid,  the  exponent  /y\.   equals  one,  and  the  equation 
becomes : 

r  =  &rt 

or 

(2.5)      y    —   Bt 

where  the  constant  B  is  the  reciprocal  of  the  Newtonian  Viscosity,  /Tl    . 
For  a  real  material  the  exponent /TV.  can  vary  between  zero  and  one,  and 
is  a  measure  of  the  relative  importance  of  the  elastic  and  viscous  contri- 
butions to  the  mechanical  behavior.   For  values  of  the  exponent  />U  be- 
tween zero  and  one,  equation  (2.3)  indicates  a  time  behavior  in  which 
stress  depends  on  both  strain  and  strain  rate.   A  material  such  as  this 
which  shows  characteristics  intermediate  between  those  of  an  ideal  solid 
or  ideal  liquid  is  called  viscoelastic. 

In  some  situations  neither  the  strain  nor  the  strain  rate  is  direct- 
ly proportional  to  stress.   An  example  of  this  complicated  relationship 
is  the  plastic  deformation  of  metals.   In  the  absence  of  this  complexity, 
and  when  both  strain  and  strain  rate  are  linear  with  stress  as  indicated 
by  equation  (2.3)  a  material  is  said  to  display  linear  viscoelastic  be- 
havior. [2] 

Linear  viscoelasticity  might  be  expected  with  materials  whose  struc- 
ture is  not  changing  during  the  time  it  is  subjected  to  an  applied  stress, 
[11]  For  such  linear  viscoelastic  behavior,  a  law  known  as  the  Boltzmann 
Superposition  Principle  applies.   This  principle  is  illustrated  in  Figure 
1,  which  relates  to  a  creep  test  where  a  material  is  subjected  to  a  con- 
stant stress,  and  the  creep  strain  changes  with  time.   For  linear  visco- 
elastic behavior,  this  strain  is  directly  proportional  to  stress: 
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(2.6)  e(t)  =  a  J(t) 

where   J  ("tj  >  the  creep  compliance,  is  a  function  representing  the  time 
effects  of  the  viscoelastic  response.   For  a  stress  CTQ     applied  at  zero 
time,  the  strain  as  a  function  of  time  is: 

(2.7)  e(t)    -    ae  3(t) 

For  some  higher  stress  (X0  +  (Tl     applied  at  zero  time,  the  strain  would  be 
increased: 

(2.8)  £(t)  -  (cr0  +  cr,)-  J(0 

Now  if  CTQ     is  applied  at  time  X.  =  O  ,  and  stress   0",  is  added  at  a 
later  time  L  =  "t|   the  strain  becomes: 


(2.9) 


e(t)  =       cr0-J(t)   +   <r,-J(t-t() 


This  is  illustrated  in  Figure  1  where  linear  superposition  of  the  strains 
at  their  elapsed  times  is  indicated.   The  stress  history  may  be  general- 
ized by: 

(2.10)       e(t)   =     IT  (Tf3(t-ti) 

tL  =  -00 

For  a  continuous  stress  history  starting  at  zero  time,  the  Boltzmann 
Superposition  Principle  can  be  expressed  in  integral  form: 

e-t 

(2.ii)       e(t)   =        Z(t-e)<r(e)de 

Je=o 
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where  Q     is  current  time  during  the  loading  history,  where  CF(Q)      is 

written  for    rl  6       '  anc*     ^s  constant  in  tne  integration.   The  creep 
strain  in  equation  (2.11)  is  therefore  represented  by  a  convolution  in- 
tegral where  (T(<9)   is  the  "forcing  function,"  and  3  f*t""©7  is  the  "re- 
sponse function. " 

In  a  stress  relaxation  test,  a  tensile  strain  is  rapidly  applied  and 
the  stress  changes  with  time.   For  linear  viscoelastic  behavior, 

(2.i2)        cr('t)   -    e-Er(t) 

Here  the  viscoelastic  time  behavior  is  given  by  E:r(.L  ;  ,  called  the  re- 
laxation modulus.   Equations  analogous  to  (2.7)  through  (2.10)  can  be 
written  for  successively  applied  strains,  so  that  for  a  continuous  strain 
history, 

e  =  t 

(2.i3)        o-(t)     =     J   £r(t-e)e(eN)  de 

6-0 

It  has  been  shown  by  Volterra  that  equations  (2.11)  and  (2.13)  are 
consequences  of  each  other.  [11,  12,  13]   These  equations  are  also  a 
statement  of  Boltzmann's  Superposition  Principle  and  a  definition  of 
linear  viscoelastic  behavior.   A  thorough  treatment  of  the  mathematical 
theory  and  implications  of  linear  viscoelasticity  can  be  found  in  refer- 
ence [14] . 
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3.    Linear  Viscoelastic  Models. 

Mechanical  models  are  not  necessary  to  describe  viscoelastic  phenom- 
ena, but  they  have  been  useful  in  describing  a  combination  of  Hookean  and 
Newtonian  characteristics  of  materials.   The  elastic  mode  of  stress  re- 
sponse may  be  represented  by  a  spring  (Figure  2) ,  presumed  to  behave  in 
an  ideal  Hookean  manner,  i.e.,  strain  proportional  to  applied  stress  and 
independent  of  strain  rate.   The  flow  mechanism  is  indicated  by  a  dashpot, 
(Figure  3)  assumed  to  be  filled  with  a  Newtonian  fluid  so  that  its  de- 
formation is  proportional  to  the  time  and  to  applied  stress,  and  inversely 
proportional  to  the  viscosity  of  the  fluid,  /(  .   The  behavior  of  a 
material  displaying  flow  superimposed  upon  elasticity  can  be  represented 
by  a  combination  of  a  dashpot  and  a  spring.   The  two  basic  models  are  the 
Maxwell  Element  (Figure  4)  and  the  Voigt  Element  (Figure  5) .   The  response 
of  these  units  to  an  applied  stress  is  illustrated  in  Figure  6.   In  order 
to  represent  a  real  material  an  infinite  number  of  these  elements  may  be 
required.   Generalized  Voigt  and  Maxwell  models  are  shown  in  Figures  7  and 
8.   Poincare  [16]  indicates  that  a  phenomenon  representable  by  any  one 
model  is  also  representable  by  an  infinite  number  of  other  models.   If 
one  dashpot  in  the  generalized  Maxwell  model  has  infinite  viscosity 
(  syj    =    co     )  ancj  if  one  dashpot  in  the  generalized  Voigt  model  has  zero 
viscosity  (  71    =  O    )  then  the  models  become  equivalent  in  their  representa- 
tion of  viscoelastic  behavior. 

The  creep  behavior  of  most  materials  can  be  considered  as  the  super- 
position of  three  stress  responses.   (Figure  9)  The  instantaneous  elastic 
strain  can  be  represented  by  a  Hookean  spring,  and  the  non-recoverable 
viscous  flow  by  a  Newtonian  dashpot.   The  Voigt  Element  is  particularly 
useful  in  describing  a  delayed  elastic  response.   The  viscous  element 
coupled  in  parallel  with  the  elastic  element  represents  a  damping  resistance 
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Figure  5 
Voigt  Element 
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Figure  6 
Response  of  Viscoelastic  Models  to  Stress 

(a)  Hookean  Solid 

(b)  Newtonian  Liquid 

(c)  Maxwell  Element 

(d)  Voigt  Element 
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Figure  9 
Creep  Behavior  of  a  Viscoelastic  Material 
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The  Four-Element  Model  For  Creep  Behavior 
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to  the  establishment  of  elastic  equilibrium.  [10]   A  corresponding  mech- 
anism in  a  polymeric  solid  might  be  the  uncoiling  of  long  polymer  chains 
from  their  minimum  energy  unstressed  positions,  which  involves  diffusion 
of  chain  segments  and  takes  considerable  time.   From  the  equations  re- 
presenting the  components  of  the  Voigt  element, 


The  solution  to  this  governing  differential  equation  gives  the  strain- 
time  relation: 


-i-.t 


(3.2)       Ed.  -    -£  ( i  -jeJn 


The  creep  behavior  of  linear  viscoelastic  materials  illustrated  in 
Figure  9  may  then  be  described  by  the  four-element  model  of  Figure  10. 
The  resulting  equation  becomes: 


v  ♦£(.'-*-*)♦£« 


(3.3)        e    =  ,—- 


where  T  -  ^>/Ei|  *s  fcne  retardation  time,  the  time  required  for 
the  delayed  strain  to  reach  a  value  )  —  /j^  ,  or  63.2%  of  its  final 
value. 

The  four-element  model  has  been  used  successfully  by  Gearhart  and 
Kennedy  [17]  in  their  study  of  cellulose  acetate  butyrate  plastics,  and 
by  Dienes  [18]  to  correlate  data  on  creep  of  polyethylene  at  elevated 
temperatures.   Marin  and  Pao  [19]  have  extended  the  four  parameter  equa- 
tion to  account  for  non-linear  stress  behavior  frequently  encountered  in 
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the  delayed  strain  and  viscous  flow  components  for  a  simple  plastic. 
Their  equation  (in  its  simplest  form)  can  be  written: 

<3.4)        e  =  |  +  c,tf~'(i-.£Clt)  +  c^ffmt-t 

A  modified  version  of  this  equation  was  developed  to  characterize  the 
creep  behavior  of  the  materials  studied  in  this  investigation. 
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4.   Non-Linear  Viscoelastic  Behavior. 

In  a  previous  study  [6]  involving  polymeric  materials  in  the  same 
loading  ranges  used  in  this  investigation,  marked  non-linear  behavior  was 
encountered  in  the  elastic,  delayed  elastic,  and  flow  regions  of  stress 
response.   Such  deviations  from  linear  viscoelasticity  are  not  uncommon. 
[20]   In  a  study  of  creep  in  cellulose  acetate,  Haward  [21]  found  that  the 
delayed  elastic  strain  component  increased  by  a  factor  of  six  when  the 
stress  was  increased  30%.   This  non-linear  behavior  occurred  at  a  compara- 
tively low  stress  of  500  psi.   Similar  stress  non-linearities  have  been 
observed  for  glass-like  amorphous  polymers  [22],  polymers  in  the  transi- 
tion region  just  below  T   [2],  and  polymers  in  the  rubbery  state  [23]. 

O 

In  all  cases  where  deviations  from  linearity  occur,  an  increase  in  stress 
produces  a  more-than-proportionate  increase  in  strain. 

Alfrey  [10]  demonstrated  that  for  large  values  of  strain,  non-linear- 
ity can  be  expected.   From  the  theory  of  elasticity  the  work  required  to 
distort  a  solid  to  a  given  strain,  called  the  "free  energy  of  deformation,1' 
can  be  represented  by  a  power  series: 


(4.1) 


5T(Y)     =     cj^  t  c,Y3  ♦  c4Y' 


The  shear  modulus  Cx     has  a  magnitude  given  by: 

(4.2)      G   =   ^Tyz 

For  small  values  of  strain  only  the  quadratic  term  is  significant  and 
the  modulus  is  a  constant,  independent  of  strain   fl   .   At  larger  strains, 
when  higher  order  terms  become  important,  the  shear  modulus  becomes  a  func- 
tion of   Y"    ,  and  stress  non-linearities  result.   A  similar  analysis 
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predicts  non-Newtonian  flow  of  viscoelastic  materials  at  higher  stresses. 

For  glassy  polymers  which  have  high  moduli  and  can  support  large 
stresses,  non-linear  behavior  has  been  observed  even  when  the  strains  are 
relatively  small.  [2]   This  suggests  that  internal  structural  changes  are 
responsible  for  the  deviations.   For  example,  it  is  reasonable  to  suppose 
that  the  flow  rate  will  be  markedly  different  for  polymer  chains  which 
are  straight  than  for  a  coiled  chain.   Experimental  measurements  in  this 
field  are  sparse,  but  some  theoretical  aspects  have  been  discussed  by 
Schwarzl  in  terms  of  transition-state  kinetics.  [24] 

Nonlinearity  of  viscous  and  elastic  elements  of  a  model  create 
systems  of  nonlinear  differential  equations  which  are  difficult  to  solve 
by  ordinary  methods.   In  this  investigation  an  equation  similar  to  that 
of  Marin  and  Pao  (Equation  3.4)  was  used  to  fit  creep  data  for  various 
polymers.   By  making  some  simplifying  assumptions  in  the  equation,  and 
with  the  aid  of  a  digital  computer,  satisfactory  predictions  of  composite 
behavior  were  obtained. 
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5.    The  Experiment. 

The  viscoelastic  materials  selected  for  the  investigation  were: 
high  density  linear  polyethylene,  isotactic  polypropylene,  nylon  6/10, 
and  polycarbonate.   All  were  obtained  from  commercial  sources  in  the  form 
of  extruded  sheet.   General  mechanical  properties  of  these  plastics  deter- 
mined by  tension  testing,  along  with  approximate  values  of  other  physical 
constants  from  the  literature  [15,  25]  are  tabulated  in  Appendix  A. 
Creep  test  specimens  were  cut  and  machined  to  the  dimensions  shown  in 
Figure  11.   The  thickness  of  the  specimens  was  that  of  the  sheet  materi- 
al; approximately  0.060  inches  for  polycarbonate  and  nylon,  and  0.192  inches 
for  polypropylene  and  polyethylene.   Specimens  used  for  relaxation  tests 
and  tensile  tests  were  made  with  the  dimensions  shown  in  Figure  12. 
All  specimens  were  cut  with  the  long  axis  parallel  to  the  extrusion  direc- 
tion. 

Creep  tests  were  carried  out  on  constant  load  creep  machines  in  a 
laboratory  with  the  temperature  maintained  at  23  +  1°  Celsius  and  relative 
humidity  controlled  at  50  +  2%.   Construction  features  of  the  machines 
and  the  method  of  gripping  the  specimen  are  discussed  in  Appendix  G. 
Series  of  creep  tests  were  run  at  different  stress  levels,  ranging  be- 
tween 10  and  80  percent  of  the  nominal  yield  strength.   Extension  of  the 
test  specimens  over  a  two-inch  gauge  length  was  measured  with  a  Tripoli- 
tis  type  extensometer  which  could  be  read  to  +  0.0001  inches.   After  load- 
ing the  specimens,  readings  were  taken  at  elapsed  times  of  6,  12,  30,  60, 
90,  120  seconds,  and  then  in  increasing  intervals  from  1  to  10  minutes 
during  the  first  hour,  then  every  30  minutes  during  the  first  day.   The 
test  was  continued  for  100  to  150  hours,  with  an  average  of  50  data 
points  taken  during  the  run.   Similar  test  procedures  were  followed  for 
the  composites,  with  the  component  materials  clamped  together  outside  the 
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Creep  Test  Specimen 
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Figure  12 
Tensile  Test  Specimen 


Figure  13 
Clamped  Composite 
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gauge  length  to  insure  uniform  elongation  between  the  clamps.  (See  Figure 
13) 

Uniaxial  tension  tests  were  run  at  the  same  temperature  and  humidity 
on  an  Instron  testing  machine.   Specimens  were  pulled  at  a  constant  rate 
of  deformation,  and  Tripolitis  extensometers  were  used  to  measure  the 
extension  over  a  two-inch  gauge  length.   Preliminary  experiments  measur- 
ing the  strain  at  timed  intervals  indicated  that  the  strain  rate  was  es- 
sentially constant  throughout  the  test. 

Stress  relaxation  tests  were  performed  on  the  Instron  tester  in  the 
following  manner.   Samples  were  strained  rapidly  at  ten  inches  per  minute 
to  a  predetermined  load  level,  where  a  load  limit  microswitch  automati- 
cally stopped  crosshead  motion.   Graphical  output  of  load  relaxation 
versus  time  at  constant  strain  was  obtained  from  the  recorder  installed 
on  the  testing  machine.   Initial  straining  was  essentially  instantaneous; 
less  than  one  second  was  required  to  reach  the  desired  stress.   Dura- 
tion of  the  relaxation  tests  was  approximately  thirty  hours. 

A  portion  of  the  creep  data  for  polypropylene  and  polyethylene  was 
obtained  in  a  previous  study.   Because  this  work  [6]  was  done  in  the 
same  laboratory  utilizing  identical  test  procedures  and  specimens  cut 
from  the  same  stock,  no  distinction  is  made  in  the  data  tabulation  in 
Appendix  A. 
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6.    Selection  of  an  Equation  to  Characterize  Creep  Behavior. 

The  four-element  model  (Figure  10)  discussed  in  Section  3  describes 
the  creep  process  in  three  stages;  an  instantaneous  elastic  strain,  a 
delayed  elastic  response,  and  non- recoverable  plastic  flow.   While  there 
is  considerable  theoretical  justification  for  this  description  of  creep 
behavior  [11,  27],  this  simple  form  of  the  model  has  two  serious  short- 
comings.  First,  it  fails  to  account  for  the  non-linear  phenomena  fre- 
quently encountered  in  tests  made  at  stresses  above  25%  of  the  yield 
strength.   Improved  fits  with  actual  data  can  be  achieved  by  modifying 
the  basic  equation  (Equation  (3.3))  to  account  for  these  stress  anomalies. 
Such  a  modification  is  shown  in  equation  (3.4),  and  this  method  was  uti- 
lized for  data  fitting  in  this  investigation. 

«•»    e-i*£(.i-i^t*;£-t 

The  second  limitation  imposed  by  the  four-element  model  is  more  funda- 
mental.  The  presence  of  only  one  Voigt  Element  implies  that  the  delayed 
elastic  response  of  a  polymer  is  a  single  molecular  mechanism  with  one 
characteristic  retardation  time.  As  pointed  out  by  Alfrey  [15],  the  coils 
and  linkages  of  a  polymer  chain  have  a  multitude  of  sizes  and  spatial  ar- 
rangements.  One  would  expect  local  kinks  to  be  straightened  out  more 
rapidly  than  long  range  coils,  and  that  elastic  equilibrium  would  be 
achieved  sooner  if  the  chain  segment  had  a  favorable  orientation  with  the 
stress  axis.   The  delayed  elastic  response  of  a  viscoelastic  material 
would  then  consist  of  several  mechanisms  with  different  retardation  times. 
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The  "memory111  effects  observed  in  high  polymers  are  frequently  cited  as 
proof  of  the  existence  of  more  than  one  retardation  time.   Failure  of  a 
single  Voigt  Element  to  represent  accurately  the  delayed  elastic  response 
in  the  materials  used  in  this  study  is  discussed  in  Appendix  F. 

Figure  14  illustrates  schematically  the  shape  of  a  creep  curve  for 
a  viscoelastic  substance  compared  to  a  set  of  simple  exponential  curves 
of  the  form  €"  =  K.  (  \  —  JL  )      ,  and  shows  the  failure  of  a  single 

time  constant  concept  to  represent  actual  data.   Figure  15  shows  the  re- 
tardation spectra  corresponding  to  the  simple  exponentials.   In  this  figure, 
the  term  }d      represents  the  rate  constant  (=  /T    )  and   Oryy\j)    the  rate 
constant  distribution  function.   The  creep  curve  of  Figure  14  can  be  ap- 
proximated very  closely  by  a  finite  number  of  simple  exponentials  appropri- 
ately weighted.   The  model  representing  the  delayed  elastic  response  would 
then  be  the  generalized  Voigt  Model  of  Figure  7.   Such  a  representation 
consisting  of  exponentials  Ci-'  ,  Q)   ,  and  v3)  is  given  in  Figure  16. 
A  generalized  equation  for  /H^  discreet  relaxation  times,  with  the  in- 
stantaneous strain  (A.)  and  the  flow  strain  (A_  x  t)  included  would  be: 

(6.i)       e  =  At   +  Aft    4    ^'(l-fat^*    ) 

where  OL'L  —   ^(bi.)  ,  the  discontinuous  distribution  function  for  the  l~" 
exponential,  and  cL^i=   1.   Although  this  "spike  distribution"  model 
was  not  used  in  this  study  to  correlate  creep  data,  G.  F.  Kinney,  in  a 
parallel  investigation  at  NPGS,  has  used  the  model  to  fit  the  same  data. 
In  his  work,  rate  constants  were  arbitrarily  assigned,  then  appropriately 
weighted  using  a  least  squares  curve  fitting  technique  adapted  to  a  digi- 
tal computer.   Two  to  six  retardation  times  were  used,  and  excellent  fits 
were  obtained,  even  with  the  simple  two-spike  distribution. 
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Figure  14 
Comparison  Of  Creep  Curve  With  Simple  Exponential  Curves 
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Figure  15 

Retardation  Spectra  Corresponding  to 
the  Simple  Exponential  Curves  of  Fig.  14 


35 


£d 


waiqhW  sum  of  exponan1'\ois  ©,(D,® ,  R9.I4 
actual  crcap  curve 


t\mc 


#) 


b, 


bt 


J*-/* 


Figure  16 

Representation  of  Creep  Behavior  By  Spike 
Distribution  of  Retardation  Times 
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Another  method  of  representing  delayed  elastic  behavior  is  by  a 
continuous  distribution  of  retardation  times,   a (  fo  ) .   This  concept  is 
attractive  from  a  theoretical  standpoint  if  one  considers  delayed  elasti- 
city as  consisting  of  many  molecular  mechanisms  whose  rate  constants  form 
a  quasi-continuous  spectrum.   Alfrey  [10]  develops  a  molecular  theory  of 
delayed  response  and  arrives  at  a  theoretical  distribution  of  retardation 
times  for  a  monodisperse  amorphous  linear  polymer.   This  distribution, 
whose  upper  and  lower  limits  can  be  determined  from  the  theory,  is  shown 
in  Figure  17.   (The  notation  here  is  left  in  the  form  given  by  Alfrey) 
For  a  polymer  consisting  of  two  different  molecular  weights,  Alfrey  has 
deduced  the  retardation  spectrum  shown  in  Figure  18.   Tobolsky,  Andrews, 
and  Dunnel  [28]  used  a  box  distribution  (Figure  19),  to  correlate  stress 
relaxation  data.   The  limits  of  such  a  box  distribution  may  be  obtained 
graphically  as  shown  in  Figure  19  from  a  data  plot  of  delayed  elastic 
strain  versus  logarithm  of  time. 

The  calculation  of  the  retardation  spectrum  for  an  actual  material, 
where  the  time  constants  extend  over  many  decades  of  time,  can  be  an  in- 
volved process.   A  simple  method  [15]  of  approximating  the  spectrum  is 
given  by: 


(6.2) 


L(r)  =  ferOM-W 


Here,[  3  (x.)  —   '?l  i  corresponds  to  the  delayed  strain  response, 
£<*  (^t )   ,  and   L  (T)   is  the  distribution  of  retardation  times  which 
can  be  converted  to  the  rate  constant  distribution  <l(b).   Figure  20 
demonstrates  this  graphically  by  plotting  the  slope  of  the  delayed  strain 
vs.  logtime  curve  versus  logtime.   It  is  apparent  that  if  Figure  20  does 
represent  the  response  of  the  four-element  model  of  Figure  10,  the  approxi- 
mation gives  a  curve  with  a  peak  at  the  most  probable  ^instead  of  a  single 
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Retardation  Spectrum  for  a  Monodisperse 
Amorphous  Linear  Polymer  (Ref.  [lo]  ) 


m 


Poly  mar  A  (lowm.w) 


L(A) 


Polymer  B  Chi^Vi  m.w.) 


Figure  18 

Spectrum  for  a  Polymer  Consisting  of 
Two  Different  Molecular  Weight  Fractions 
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Figure  19 
Box  Distribution  Of  Retardation  Times 


Figure  20 
Method  Of  Approximating  Retardation  Spectra 
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spike.   For  real  materials,  however,  the  approximation  is  better.   Refer- 
ence [27]  reviews  more  accurate  methods  for  estimating  the  retardation 
spectra. 

Analysis  of  the  box  distribution  function  is  given  in  Appendix  F. 
It  is  shown  that  this  distribution  is  unsuitable  for  characterizing  the 
creep  behavior  of  the  materials  tested,  as  are  other  monomodal  distribu- 
tion functions. 

In  the  absence  of  a  mathematically  tractable  function  to  represent 
the  retardation  spectrum  of  the  materials  tested,  the  continuous  spectrum 
concept  was  discarded.   Instead,  a  simple  empirical  equation  was  developed 
by  assuming  that  the  rate  constant  associated  with  the  delayed  elastic 
response  is  a  function  of  time.   One  might  expect  that  a  polymer  chain 
segment  should  exhibit  a  changing  retardation  time  as  the  uncoiling  pro- 
cess continues  and  the  shape  of  the  molecule  changes.   Thus,  an  overall 
characteristic  time  constant  7r   might  be  considered,  with  a  simple 
power  dependency  on  time : 

(6.3)        TR  »  CLi 

Substitution  of  this  expression  into  the  equation  for  the  four  element 
model  (Equation  (6.4))  leads  to  an  equation  of  the  form  of  equation  (6.5) 

(6.4)  g  =  A*    +  Ait     +    A3l\-^~t/rR) 

(6.5)  e=  A,   +■  Ait     +    A3(\-X        ) 
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Equation  (6.5)  was  used  to  fit  creep  data  in  this  investigation.   The 
method  of  obtaining  the  creep  parameters  A..  ,  A-,  A_,  C,  and  Q  are  dis- 


cussed in  the  next  section. 
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7.   Analysis  of  the  Data. 

Creep  data  were  taken  for  all  specimens  in  the  form  of  elongation 
versus  time.   These  data  were  punched  on  cards  and  a  print-out  obtained 
from  an  IBM  360  Digital  Computer.   The  printout  served  as  a  check  to  pick 
up  gross  errors  due  to  equipment  malfunction  during  the  testing  period  or 
errors  made  in  transcription  of  the  data. 

The  first  step  in  fitting  the  creep  data  was  to  obtain  the  time- 
independent  part  of  the  elastic  strain.   This  instantaneous  elastic 
strain,   £-o  ,  occurring  when  the  load  is  applied,  is  represented  by  the 
parameter   Aj  in  equation  (7.1): 

(7.D         e  =  A,    +  Att  +  Ai(l--e"      ) 

To  evaluate  the  constant  f\\    ,   a  direct  method  was  employed.   The  first 
four  data  points  were  used  to  perform  an  extrapolation  back  to  zero  time. 
This  extrapolation  to  determine  the  parameter  /\\    is  illustrated  in 
Figure  21.   It  is  clear  from  this  figure  that  this  evaluation  of  Al 
is  at  best  a  rough  estimate  due  to  the  nonlinear  nature  of  the  extrapola- 
tion.  It  will  be  shown  in  Appendix  D  that  the  parameter  A|   does  in- 
clude a  portion  of  the  delayed  elastic  strain,  but  for  the  time  frame  of 
this  experiment,  the  strain   r\j  occurs  essentially  instantaneously. 

To  account  for  the  non-linear  viscoelastic  behavior  encountered,  the 
parameters   A(  ,   A^  ,  and  A3  were  correlated  with  stress  so  that 
the  creep  behavior  could  be  predicted  for  any  stress  level.   The  general- 
ized form  of  equation  (7.1)  becomes: 


<7-2>       e=  KcrN    +    baRt    +   PffM(i-Act*) 
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Figure  21 
Extrapolation  of  Creep  Data  to  Determine  Instantaneous  Strain 
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The  magnitude  of  parameters  A7.  and  A3  depend  on  the  value  chosen  for 
/\  (   .   Prior  to  evaluating  A 1.    and  A3    ,  the  extrapolated  values 
of  zero  time  strain   A|  were  correlated  with  the  first  term  of  equation 
(7.2): 


(7.3)         A,    =    KtfM 


(7.4)     J^A\     -      JUK      +■  N  1UC 

Values  of  Jl/v^  A\\     for  different  stress  levels  were  read  into  a  digital 
computer,  and  equation  (7.4)  was  used  to  perform  a  least  squares  fit  to 
determine  the  stress  coefficient   rx   and  the  stress  exponent  N 
Figure  22  shows  the  plot  of  equation  (7.4)  and  the  fitted  curve.   The 
data  displayed  in  this  figure  are  that  obtained  from  the  extrapolations 
of  Figure  21,  made  for  polycarbonate. 

An  independent  method  of  obtaining  the  zero-time  strain  was  attempted 
by  running  a  series  of  tensile  tests  at  various  strain  rates,  and  plotting 
the  resulting  data  as  sketched  in  Figure  23.   From  these  stress-strain 
curves,  a  cross-plot  of  strain  versus  reciprocal  strain  rate  may  be  con- 
structed as  illustrated  in  Figure  24.   The  zero- time  strain  can  then  be 
obtained  by  extrapolating  the  data  to  strain  rate  £  =  &&        ,  which 
corresponds  to  instantaneous  loading.   This  technique  was  attempted  using 
a  digital  computer  to  perform  a  parabolic  extrapolation,  but  the  exponential 
character  of  the  data  plot  invariably  resulted  in  high  estimates  of  the 
zero- time  strain,  So      ,  as  indicated  in  Figure  24.   In  addition,  the 
maximum  deformation  rate  available  on  the  Instron  tester  used  for  the 
tensile  tests  was   £  =  10  inches /minute.   Consequently,  there  was  a  lack 
of  data  in  the  region  of  greatest  interest,  i.e.,  where  /l  — **  0.   This 
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Figure  22 
Correlation  of  Parameter  A^  with  Stress 
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Figure  23 
Stress-Strain  Curves  for  Various  Strain  Rates 
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Figure  24 
Zero  Time  Strain  Extrapolations  from  Stress-Strain  Data 
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method  may  be  expected  to  produce  excellent  results  if  high-speed  testing 
apparatus  were  available  [29],  but  in  view  of  the  limitations  discussed 
above  it  was  discarded  in  favor  of  the  creep  data  extrapolation  technique. 

After  the  parameter  A.,  was  found  by  extrapolation  and  "smoothed"  by 
the  stress  correlation  (Equation  (7.3))  the  creep  data,  along  with  the 
calculated  parameter  A- ,  was  placed  in  the  digital  computer  via  program 
FITDATA  (See  Appendix  B) .   The  viscous  flow  parameter  A„  was  determined 
by  a  linear  extrapolation  using  data  points  at  elapsed  times  greater  than 
70  hours.   It  was  assumed  that  elastic  equilibrium  was  achieved  by  this 
time.   With  parameter  A.  specified  and  parameter  A-  computed,  the  maximum 
magnitude  of  the  delayed  elastic  strain,  A^ ,  is  determined.   This  pro- 
cedure  used  by  program  FITDATA  is  outlined  in  the  sketch  of  Figure  25. 
The  time  constant  parameters  were  computed  by  FITDATA  by  evaluating  the 
delayed  strain  fraction,   £.  ,  at  various  times. 

(7.5)    £  =  e*/A,  =  I  -  4-"ctQ 

from  which  ,  ^ 

-jMi-£)   =  ct 

(7.6)       /*,[-,Mt-&)]  =  z*-c  +  o.u,t 

A  plot  of  the  left-hand  side  of  equation  (7.6)  versus  XA^L.       was  given 
a  least  squares  curve  fit,  and  the  constants  C  and  Q  determined.   The 
method  is  sketched  in  Figure  26,  and  graphical  computer  output  of  actual 
data  is  displayed  in  Appendix  B. 

The  creep  parameters  A  and  A«  were  correlated  with  stress  using 
equations  (7.7)  and  (7.8): 

a.7)        At    =     5crR 
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h*  t>TOhrs. 


Figure  25 
Determination  of   Parameters  k%  and  A* 


ln[-b(i-£)] 


Figure  26 
Determination  of  Parameters  C  and  Q 
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(7.8)      A3   =   P^^ 


The  method  used  for  the  A.  correlation  was  employed  to  determine  the  stress 
exponents  and  coefficients  in  equations  (7.7)  and  (7.8).   These  determina- 
tions for  polycarbonate  are  illustrated  in  Figures  27  and  28. 

A  question  arises  as  to  the  validity  of  inclusion  of  the  viscous 
flow  term  A2/L  in  equation  (7.1).   If  the  term  is  not  included  and 
the  first  and  second  derivatives  of  strain  with  respect  to  log  time  are 
taken,  at  long  times  the  second  derivative  is  always  negative  and  becomes 
vanishingly  small-   At  long  times  the  first  derivative,  corresponding  to 
the  slope  of  the  6.  vs.  J^y\^-     curve  of  Figure  29,  approaches  zero.   With 
no  viscous  flow,  elastic  equilibrium  is  eventually  reached  and  extension 
ceases.   If  the  flow  term  A^*  L    is  included,  the  second  derivative 
changes  sign  and  becomes  increasingly  positive  at  long  times.   This  behav- 
ior is  observed  as  an  inflection  and  an  upswing  in  the  plot  of  £    vs.  j2W"t-  , 
Figure  30.   If  viscous  flow  is  a  factor  in  creep  behavior,  therefore,  a 
data  plot  of  £  vs.  XM,   u   should  show  the  inflection  and  upturn  indi- 
cated by  Figure  30.   In  all  the  materials  investigated  this  behavior  due 
to  viscous  flow  was  observed,  as  illustrated  by  the  example  of  Figure  31, 
justifying  inclusion  of  the  term   A  2.  *  "t    in  equation  (7.1). 

Values  of  the  creep  parameters  determined  for  the  materials  tested 
are  summarized  in  Table  2  of  Appendix  A.   The  constants  C  and  Q  are  ap- 
parently independent  of  stress  and  are  of  the  same  magnitude  for  all 
materials  tested.   Standard  deviations  in  strain  were  acceptable  for  all 
runs.   The  errors  in  computed  strain  were  less  than  5%  in  the  region  where 
delayed  elastic  behavior  is  significant,  except  at  elapsed  times  under 
two  minutes.   Sample  computer  output  from  F1TDATA  listing  point  by  point 
comparison  of  observed  and  computed  strain  is  given  in  Appendix  B.   For 
typical  values  of  the  constants   C  =  0.8  „  and   Q.  =  0."35  ,  the  delayed 
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Figure  28 
Correlation  of  Parameter  Ag  with  Stress 
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Figure  29 
Strain  versus  Ln  Time  with  No  Viscous  Flow 


Figure  30 
Strain  versus  Ln  Time  with  Viscous  Flow 
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strain  is  only  97%  exhausted  at  an  elapsed  time  of  70  hours.   For  the 
extrapolation  to  determine  parameter  A9,  elastic  equilibrium  is  assumed. 
Consequently,  in  the  resulting  "forced  fit,"  after  the  other  parameters 
are  determined,  the  computed  strain  is  approximately  one  to  two  percent 
low  at  long  times. 

From  the  tests  at  higher  stresses,  it  is  apparent  that  elastic  equi- 
librium is  not  achieved  during  the  testing  period.   This  produces  a  high 
estimate  of  the  parameter  A_,  which  in  turn  results  in  a  low  estimate  for 
parameter  A~ .   At  low  stresses  the  flow  rates  are  very  small  and  the  er- 
ror introduced  is  negligible.   However,  power  dependence  of  A_  on  stress 
creates  noticeable  deviations  at  high  stress  levels.   This  deviation  is 
clear  in  Figures  27  and  28  at  stress  levels  of  CX  =  4500  psi  and  (X  = 
5000  psi.   In  performing  the  stress  correlations,  data  at  lower  stresses 
was  weighted  more  heavily,  which  corrects  the  erroneous  values  of  A  and 
A3-   (See  figures  27  and  28). 

The  fitted  empirical  equation  (7.2)  provides  excellent  predictions 
of  creep  behavior  on  individual  materials  loaded  at  constant  stress. 
When  a  simple  composite  sandwich  of  two  materials  is  tested  in  creep  how- 
ever, the  stress  in  each  component  changes  with  time  due  to  the  differ- 
ent viscoelastic  responses  of  each  material.   If  the  time  responses  were 
linear  in  stress,  accurate  theoretical  predictions  of  composite  behavior 
could  be  made.   Reference  [15]  outlines  procedures  for  calculating  stress 
as  a  function  of  strain  and  time  when  the  creep  compliance  j("t)  is 
linear  in  stress.   The  Boltzmann  Superposition  Principle  (Equation  (2.11)) 
could  then  be  applied  to  determine  the  creep  behavior  of  a  composite. 
From  Table  2  in  Appendix  A,  the  stress  exponents  are  greater  than  one,  so 
the  viscoelastic  behavior  of  the  materials  tested  is  non-linear  and  the 
creep  compliance  J  ft)  *s  not  independent  of  stress.   Resulting  stress- 
strain-time  relationships  become  extremely  complex.   If  a  non-linear 
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viscoelastic  "equation  of  state"  could  be  developed,  it  might  be  too  compli- 
cated for  practical  use,  even  with  digital  computer  techniques. 

In  order  to  predict  the  behavior  of  composites  studied  in  this  investi- 
gation, a  simplifying  assumption  was  made.   The  strain  response  of  a 
material  was  assumed  to  be  independent  of  its  stress-time  history.   Equa- 
tion (7.2)  then  functions  as  an  equation  of  state;  if  a  stress  and  time 
are  specified,  then  the  strain  is  determined.   From  the  physical  restric- 
tions placed  on  the  composite  tested  in  creep, 

(7.9)  P^  +     Pb   «P 

(7.10)  Get  =  €b  ~  e 

where  P  represents  applied  load  and  the  subscripts  a  and  b  are  used  to 
identify  quantities  pertaining  to  component  A  or  component  B  of  the  com- 
posite.  Letting  a  and  b  denote  the  relative  fraction  of  component  A 
and  component  B  respectively, 

(7.11)  OL  *  b   =■   I 

(7.12)  a  01  +•   bcrb   ==  <T 

Unsubscripted  symbols  above  pertain  to  the  composite  structure  as  a  whole. 
From  equation  (7.2), 

(7.13)       €b=  KbffbNb+  BkCt  t  P,crbMb(\-/Cbt  ) 
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Equation  (7.12)  may  be  solved  for   (Xl  in  terms  of  known  composite  stress 

(J        ,  known  volume  fractions   a  and  b,  and  the  unknown  stress  0"^  . 
With  this  result  substituted  for  Go.    in  equation  (7.14),  equations  (7.13) 
and  (7.14)  may  be  equated  according  to  equation  (7.10).  For  a  specified 
time,  the  corresponding  value  of  £Tb  may  be  found.   A  series  of  values 
of  (J^    and  time  can  then  be  substituted  in  equation  (7.13)  to  determine 
strain  £^  according  to  the  equation  of  state  assumption.   By  equation 
(7.10),  strain   G^  is  the  same  as  the  strain  in  the  composite. 

A  digital  computer  program,  program  COMCUR,  was  used  to  predict  compos- 
ite behavior  in  this  manner.   Sample  output  from  COMCUR  consisting  of  point- 
by-point  comparisons  of  predicted  and  observed  strain  is  found  in  Appendix 
C,  along  with  graphical  displays  of  this  output.   A  summary  of  results  of 
this  prediction  method  are  tabulated  in  Table  4  of  Appendix  A.   When  the 
viscoelastic  responses  are  similar,  the  load  distribution  in  the  components 
of  the  composite  does  not  undergo  large  changes  during  the  test  and  pre- 
dictions are  excellent.   Such  is  the  case  for  the  polyethylene/polypro- 
pylene (PE/PP)  composites.   It  should  be  noted  that  if  the  stress  in  the 
components  were  constant,  the  method  would  be  theoretically  correct,  since 
equations  (7.13)  and  (7.14)  are  valid  for  a  constant  stress  test.   Large 
errors  in  predictions  occur  with  component  materials  of  nylon  and  poly- 
carbonate.  As  indicated  by  the  magnitude  of  the  creep  parameters  given 
in  Table  2,  Appendix  A,  the  viscoelastic  response  of  these  materials  is 
quite  different.   Graphical  output  from  COMCUR  shows  that  the  load  dis- 
tribution for  the  NY/PC  composites  changes  rapidly  during  the  first  few 
minutes  of  the  test. 

The  empirical  creep  equation  (Equation  7.2)  was  used  to  predict  stress- 
strain  behavior  and  stress  relaxation  behavior  for  single  materials.   This 
analysis  was  done  with  a  digital  computer  using  the  equation  of  state 
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assumption.   An  explanation  of  program  RELAX,  used  to  predict  stress  re- 
laxation behavior  is  given  in  Appendix  D.   Output  from  program  SIGEPS,  an 
estimate  of  stress-strain  curves,  is  given  in  Appendix  E.   In  both  cases, 
computed  responses  were  considerably  different  from  data  obtained  by  experi- 
ment. 
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8.    Conclusions. 

The  results  of  computer  programs  RELAX  and  SIGEPS  demonstrate  the 
complex  nature  of  viscoelastic  behavior.   In  tensile  tests  and  relaxation 
experiments  the  stress  changes  with  time  and  the  constant  stress  equation 
gives  poor  predictions  of  responses.  While  the  equation  gives  excellent 
estimates  of  creep  behavior,  the  stress-time  history  cannot  be  ignored 
when  predicting  viscoelastic  behavior  in  other  modes  of  testing. 

The  combined  mathematical  and  experimental  approach  to  understand- 
ing composite  behavior  is  a  powerful  tool.   An  analysis  of  simple  compos- 
ite structures  provides  guidelines  for  predicting  behavior  of  more  compli- 
cated systems  where  factors  such  as  component  distribution  and  interfacial 
shear  stresses  become  important.   The  non-linear  behavior  of  many  visco- 
elastic materials,  even  at  moderate  stresses,  leads  to  involved  mathemati- 
cal relationships.   A  digital  computer  is  necessary  for  analysis  of  these 
complex  systems. 

Some  opportunities  for  further  investigation  are: 

(1)  Use  of  X-ray  diffraction  techniques  to  follow  structural 
changes  during  testing  in  order  to  determine  the  relation 
of  internal  structure  change  to  observed  non-linear  visco- 
elastic effects. 

(2)  Fabrication  and  testing  of  more  realistic  composite  materi- 
als to  determine  if  the  empirical  equation  developed  here  is 
general  enough  to  characterize  the  behavior  of  these  materi- 
als in  creep. 

(3)  Use  of  other  modes  of  testing,  such  as  stress  relaxation  and 
dynamic  mechanical  tests  to  assist  in  development  of  non- 
linear viscoelastic  theory. 
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APPENDIX  A 
TABULATION  OF  RESULTS 

Table  1 
Properties  of  Tested  Polymers 


Material 

*Tensile 
Strength 
psi 

*Tensile 
Modulus  _ 
psi  x  10 

*Percent 
Elong. 

Melt.     Glass 

Point  Transition 

°C     T  ,   °C. 
!   g 

POLYETHYLENE 

4000 

220 

380 

137 

-120 

POLYPROPYLENE 

4400 

190 

210 

176 

-10 

POLYCARBONATE 

8800 

300 

60 

267 

150 

NYLON 

5900 

120 

185 

227 

40 

i 

1 

*At  strain  rate  =  0.012  in/in/min 
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APPENDIX  B 
COMPUTER  PROGRAM  FOR  DETERMINING  CREEP  PARAMETERS 
A  Fortran  IV  computer  program,  compatible  with  the  IBM  360  system, 
was  used  to  determine  the  creep  parameters  A_,  A_,   C,  and  Q  in  the 
equation: 

e  -  A,  -r-  Avt  +  A3(\-X~c    ; 

The  program,  called  FITDATA,  is  entered  with  the  predetermined  parameter 
A-  and  the  deformation- time  data  obtained  from  the  creep  tests.   Computer 
output  consists  of  the  calculated  parameters,  tabulated  values  of  observed 
and  computed  strain,  and  intermediate  data  used  in  the  parametric  determina- 
tion.  An  off-line  plotting  routine,  called  SUBROUTINE  DRAW,  is  used  to 
display  the  output  (experimental  data  and  fitted  curve)  in  graphical  form. 
This  subroutine  was  also  used  to  plot  the  data  for  determining  parameters 
C  and  Q  in  the  equation  above.   Least  squares  fitting  was  done  by  the 
library  subroutine,  SUBROUTINE  LSQPOL  available  at  the  NPGS  Computer 
Facility.   Program  FITDATA,  along  with  typical  tabulated  and  graphical 
output,  is  given  on  the  following  pages.   The  program,  discussed  briefly 
in  Section  7,  is  annotated  with  comment  cards. 
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COMPUTER    PROGRAM    "FITDATA" 


FlkIRAN     IV    ('.    LEVEL     1,     MUD    1  MAIN                                         RATE    =    68139                           0<>/32/3<)                                      P»r,n    cr01 

U.-C1  IMPLICIT    RPAL*B     (A-h),     REAL    *B     (O-ZI 

00C2  RtAL*H     I  TI  TLE  I  li  I  ,     KTITLEI12I,     LARA/BH                        / 

OC03  DIMENSION  tP! ICO) , T( 100) ,  J!  (  100  I  ,  E  PZ  I  1  00  I  ,  Y  (  100  I  ,  DFLYll 
1001,  tl<5),SB(6l,T6(5),ST(5l,CQ(5),SC<5),AQ(30,30),FO(100).TLN(100l 
2,t0Z(  100  1, TLN7 ( 100  1 ,TGG!200) ,YGK 100)  ,ERR ( 100) ,EPRF( 100) ,FFDI IOC)  , 
lEPO<100l,TO(100l,A(30,301,rT(lGOI1x(100l,F2(100),C(S),EP77(10D) 

uCCh  RF  AL     LABEL  ISO) ,IA8B,TG! 200) ,YG<206)  , YOM(200)  ,TSlton),f(100),TAI10C 
ll,tM(100l,(>IPIT,l)lFEP,DlST,TSKAL,EPSKAL,R,7ILCHMH             /, 
2PMKP/4H          PC,  PMARK!  1001,  TYl  100).  TLNY!  100),  TYSI  1001,  EPS!  100),  TLO, 
3FLAG,CZ,ZS<AL,YSCAL,TU(i00l , YU( 100 ) , FOLO,    W( 100) 

Ok.C5  REAL     AA1,mA2,AA3,CCC,Q0Q 

C 

C  FOR  EACH  COMPUTER  RUN  THE  FOLLOWING  INFORMATION  IS  READ  IN 

C  NDAT  =  NUMBER  OF  SETS  OF  OATA 

L  IX  =  NUMBER  OF  SETS  OF  DATA  TO  RE  DRAWN  ON  FACH  GRAPH. 

C  JC  =  CUUNTING  SYMBOL  TO  CHANGE  THE  POINT  PLOTTING 

C  SYMBOL  FOR  EACH  DATA  SET. 

C  LABEL  =  PUN  NUMBERS  USED  TO  LABEL  EACH  DRAWN  CREEP 

C  CURVE. 

C  1 T I TLE  =  TITLE  OF  GRAPH  OF  CREEP  CURVES 

C  KTITLF  TITLE  FOR  GRAPH  SHOWING  C  AND  0  DETERMINATION. 

C 

0CC6  KEAC  16,  NDAT, IX 

^CC7  16  FURMAT12I5) 

*;Ci.H  READ  17,  (LABELI  I  )  ,  I  =  l.NDAT) 

OCC<;  17  FJKMATI20A'.) 

CC1C  READ  1H, ( ITITLt ( 1 ) , I  =  1,12) 

vCll  IB  E0RMATI6AB) 

0^12  NSkIP  =  0 

l.i.13  READ  '0  2,  (KTITLEII),  1*1,121 

0C1«.  t>;2  FdRMAT(6ARI 

0C15  503  AIX  =     IX 

oClt  FNCAT  =  NDAT 

iC17  PAGL  =  FNDAT/AIX 

Oi^le  NPAGE  =  PAGE 

JCIS  NOO  =  NPAGF*IX 

C^tC  KC  =  C 

j:<;l  Of  ^1   K  =  1, NPAGE 

L- 

C  A  NEW  GRAPH  SHEET  IS  STARTED  EACH  TIMf  THF  INDFX  OF  THIS 

C  DU  LOUP  IS  CHANGED. 

L 

-Ct2  JC  =  1 

UC.J  -0  OC  <i2  L  =  1,  IX 

uc:<.  jc  =  jc  *i 

OC25  KL  =  KC  ♦  1 

^C26  PRINT  10 

GC27  10  FORMAT! 1H1 I 

C 

C  THt  EOILOWING  DATA  IS  READ  IN: 

C  NPTS  =  NUMBER  CF  POINTS  IN  THF  DATA  SFT. 

L  GL  =     TEST  GUAGE  IENGTH 

L  Al  =  PREDETERMINED  ZERO  TIMF  STRAIN 

C  TMAX  =  MAXIMUM  TIME   IN  HOURS  TO  WHICH  THE  EITTFD  CUCVF 

C  IS  TO  BE  DRAWN. 

C  tPSKAL  =  STRAIN  SCALE  ON  CREEP  CURVE  GRAPH  (IN/IN/IN(H| 

C  TSKAL  =  TIME  SCALE  ON  CREEP  CURVE  GRAPH  ( HOURS/ INCH) 

C  CONST  =  CONSTANT  USED  TO  CUNVERT  INPUT  CREEP  DATA 

C  TC  STRAIN 

f  FUN  =  RUN  DESIGNATOR  USED  TO  IDENTIFY  THE  DATA 

C  tP  =    ARRAY  OF  CREEP  DEFORMATION  DATA:  USUALLY  RFAD  I  r>: 

C  AS  DEFORMATION  IN  MRS 

L  TT  =  ARRAY  OF  TIMES  CORRESPONDING  TO  TH^  ARRAY  »"» 

C 

^026  «fc*C  23,  NPTS, f.L  ,A1  , TMAX, EPSKAL  ,  TSKAL, C  INST, RUN 

CCil  2)  FORMAT! I5,6F10.C,5X,A5> 

CCiC  REAP  2S,  (FPU),  I  =  l.NPTSI 

oC  31  25  FORMAT!  1<,ES.  II 

tC  Jt  REAO  26,  ( TTI  I  )  ,  1  =  1  ,NPTS) 

CC-)  2e  FORMAT! 7F1J. 51 

C 

C  CONVERT  THE  DEFORMATION  DATA  TO  STRAIN  IN  I NC H^ s/ | mch . 

C  "AKt  THE  INITIAL  VALUE  OF  STRAIN  EOUAI  To  THE  7FRC  TIMt 

C  STRAIN  AND  THE  INITIAL  TIMfc  EQUAL  7ER0.   NUMBER  OF  DATA 

C  POINTS  IS  INCREASED  BY  ONE,  BECOMES  "NPTS1". 

C 

...I-.,  trill   =  ai 

;  *.  3  5  r  i  i  i  i  =  o .  ■'; 

0C3c  !C    <>3       I     =    l.NFTS 

„C.-7  ill     *     1)     =    EFI  I  ) /(CCNST»GLI 

,.  C  3  r  T'.II»1)     =    TT([) 

^.C  J=i  \PTS1=NPTS«L 

..■•;  OC  30a   I  -  l^PTSl 

..-.I  -:pi  I )  =  epoi  1 1 

CC<.2  30B     TTIII     =     Tulll 
C 

L  THc     NEXT     SIXTEEN    CAPOS    SELECT    ELAPSED    TIMFS    GPFATEe     THaf 

C  70*    OF     THE     TOTAL     TESTING    TIMf.        |F    IFSS    THAN    FOUR    POINT1: 

L  Akt    OBTAINED,     THE    NFXT    LONGEST     TI"E     IS    ScLFLTFn,     fNO    THE 

C  PROCESS     IS    CONTINUED    UNTIL    FOUR     POINTS    ARE    DPT4INFP.        IF 

C  THF     SHORTEST     TIME     IS    LESS    THAN    SO*    OF     THF    I'lUI     RAPS-0 

C  T]"E,     TIE     CATA     IS    DEEMED    UNRELIABLE    FOR    DETERMINING 

C  PARAMETER    AJ,    AND     IS    DISCARDED.           THE     SELECTED    TIMf1:     FOR 

C  THf     FaTraPhLATICN    ARF     STIREO    IN    THE     A  n  R  \  y    «T7«    4NI,    THF 

C  CuRRESFuNOING    STRAINS    ARE    STORE)     IN    THF    1«1»    "EP7". 
C 

,t-,:  If     =     . .S«TT(NPT >1 I 

:C<,  IF     =    v..7»TT!NPTSl  ) 

LC5  KA    =    C 

CCft  30    XA    =    ka    .     L 

..<,7  IEITT(KA)-TF)3:<i,3lC,M. 

i_C°  31C    KA    =    KA-1 

CCl«  KB    =     NPTS1     -     K,. 

CCiO  IFUB    -    <.)     31C  ,  ill, 311 

C;;i  311     IF  I TT(KA»1) .LT. TM)        GO    TO     777 

l".32  GL.TC312 

;>.;3  777P^INT77o,    RUM 

;v.-,-,  77d  FORMAT!/// ,  1CX,  'DATA  FOR  RUN  ■  ,  A  5 ,  •  UNCEIIABIF  FOR  STFATY-STflTF  (R 
1FFP     t  XTKAPUL AT  I   JN1  , // ) 

01.55  GO    TO    ^d 

oC56  312     OC    313       I     =     1.K3 

;  .-.  I  T7(  1  )     =     TTIKA*!  I 
313     FF7(  I  )     =    tP(KA» I  I 

EXTrAF 
JPOL.     TF 

C  Ni    'lAL     FIT    TO    THF     SELECTED    POINTS.       fOcFFICItNTS    OF     THF 

C  FITTED     STRAIGHT    LINF    ARE    STORED    IN    THf     ARRAY    CJllfil    "«". 

C  THr     SUBROUTINE    PRINTS    OUT     INPUT    DATA    AND    THr     RF^UITC     fF 

C  THr;     FIT     IN     TABULAR    FORM.       THF    VALUES    OF     THF    CRFFP 

(.  PA"AMET-RS     A2    AND    A3    ARE    THEN    DET  =a  M  |  Nrn 

C 

..-,-,  I.ALL      L5  I.  PLl  (  K  u,  I,  0,  0,  J,  S  I  C.  MA,  T/,  CPJ  ,  W,Y,DflY,R,S1,  TO,  ST,  CO.^f,  A  0) 

,;!.,  IF!  5l(.MA.r,T.0.O->)       GO    TO    777 

.  C  f  1  A  2    =     I'  I  2  ) 

.11  2  A  3    =    M { 1 1     -    A 1 
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COMPUTF  THF  DELAYED  STRAIN  FRACTION  AS  ARRAY  "ED"  AND 
NATURAL  LOGARITHM  OF  TIME  AS  THE  ARRAY  "TLN" 

till  C  SM  iSlfl*-  ieP<t|N-T?Al  .  A2.TT.I.../A3 

0065  TLNI I)  =■  0.0 

0066  DO  322   I  =  2.NPTS1 

0067  322  TLN(I)  =  OLOGITTUll 

C  COMPUTE  THE  DATA  REOUIRED  TO  DETERMINE  CONSTANTS  C  AND  0. 

DISCARD  THE  LAST  TEN  POINTS  WHICH  MAY  GIVE  ED  VALUES 
GREATER  THAN  ONE,  DUE  TO  THE  "FORCED  FIT"  OF  PARAMFTFR  A2. 
ELIMINATE  ANY  VALUES  IN  THE  INITIAL  PART  nF  THE  RUN  WHFRF 

r  tO  MAY  HE  LESS  THAN  ZERO  BECAUSE  QF  AN  OVERESTIMATE  OF 

C  PARAMETER  Al.   ATORE  THE  RESULTS  OF  THF  CALCULATION  IN  THF 

C  ARRAYS  "EDZZ"  AND  "TLNZ". 


C 


10 
0 


THE  NEXT  EIGHTEEN  STATEMENTS  REDUCE  THE  DATA  POINTS  SO 
THAT  NO  TWO  DATA  POINTS  ARE  CLOSER  THAN  0.2  INCHES  APART 
ON  THE  GRAPH  SHEET 


0070  339  IP  «  IP  ♦  1 

0C71  IF  <<IP«LQl  -  NPTS1)  3*0,3*0,3*1 

0072  3*0  EDZUPI  «  ED(IP»LOI 

OC73  TLNZl  IP)  =■  TLNI  IP  ♦  LOI 

007*  IFI l.-EOZI IPI )  3*1,3*1,339 

0075  3*1  KPTS  =  I P  -  1 

0076  DO  325   I  =  l.KPTS 

CC77  EDZZ(I)  =  -OLOGd.  -  EDZ(D) 

0078  325  EOZZ(I)  =  ULOGI EOZZ I  I  I  I 
C 

C  USE  SUBROUTINE  LSOPOL  TO  PERFORM  A  LINEAR  FIT  TO  THE  DATA 

C  "EDZZ"  AND  "TLNZ".   PRINT  OUT  THE  DATA  AND  THE  RESULTING 

C  FIT. 
C 

0079  CALL  LSOPOL (KPTS, 1,0, 0,0, SIGMA, TLNZ, EOZZ.W.Y, DEL Y,B, SB, TO, ST, CO, 
1SCA0I 

CC80  00  205   I  =  l.KPTS 

0081  FED(I)  =  BUI  ♦  B(2)*TLNZII) 

0C82  205  ERRF(I)  =  HEDZZ(I)  -  FED  (  I )  I  »1  00.  I /EOZ  Z  (  I  I 

0083  PRINT  71,  RUN 

008*  71  FORMAT!///, 10X, 'THE  DATA  USEO  TO  FINO  "C"  AND  "0"  FOR  RUN  ',A5, 
U  IS  LISTED  BELOW: ',///. 6X, 'DATA  PT.  LN  TIME  LN (-LN ( 1 . -ED ) » 
2FITTED  VALUES    PCT . ERROR •,// > 

CC85  DO  210   I  =  l.KPTS 

OCot  II  =  I  ♦  LO 

0C87  210  PRINT  72,  I  I , TLNZ (  I  I , EDZZ (  1 1 , FEOI I  I  , ERRF ( I  I 

0088  72  FORMAT! I  12 , F 12 . *,F 1 3. * ,F 16.*, Fl 3. 2 ) 

CC8S  IFISIGMA.GT.0.1 I   GO  TO  B91 

OCSO  GO  TO  893 

CC91  891  PRINT  892,  SIGMA 

0C92  892  FORMAT!///, 10X, 'SIGMA  =  ',F5.3,'   UNRELIABLE  FIT  FHR  FINDING  C  AND 
1  0',//l 

C 

C  FROM  THE  COEFFICIENTS  OF  THE  FITTED  POLYNOMIAL,  COMPUTE 

C  THE  CREEP  PARAMETERS  C  AND  0. 

C 

0CS3  893  U  =  B(2I 

CCS*  CC  =  OEXPIBI 1) I 

C 

C  THE  NEXT  NINE  STATEMENTS  ASSIGN  200  VALUES  TO  THE  ARRAY  OF 

C  TIMES  "TGG"  IN  EQUAL  INCREMENTS  FROM  ZERO  TO  TMAX. 

C  VALUES  OF  FITTED  STRAIN  ARE  COMPUTED  IN  FUNCTION  SIIBROGRAM 

C  "EON"  USING  THE  DETERMINED  CREFP  PARAMETERS,  AND  STORED  IN 

C  THE  ARRAY  "YG". 

C 

0C95  R  =  TMAX/200. 

CC96  TGI1)  =  0.0 

CCS7  TGGU1  =  0.0 

CC98  00  **   I  =  2,200 

CC95  TGI  I )  =  TGI I-ll  ♦  R 

C1C0  **  TGG( I  I  =  TGI  I ) 

C1C1  00  *5   I  =    1,200 

0102  YGII)  =  E0N(A1,A2.A3,CC.Q,TGG( I  II 

C103  *5  YGMI I  I  =  -YGI I  I 

C 

C 

C 

C 

c 

OK*  E(l»  =  EPU) 

01C5  TAI 1)  =  TTI1I 

01C6  JA=  1 

01C7  IZ  =  1 

01C8  KR  =  2 

01C9  60  JA  =  JA  *     1 

CUC  IFI  JA.GT.NPTSI  GO  TO  62 

GUI  TSCAL  -    TSKAL 

0112  ESCAL  =  EPSKAL 

0113  D1FFT  =  (TT(JA)  -  TT ( I Z ) I /TSC AL 
CU*  DIFEP  =  (EP(JA)  -  EPI  IZII/ESCAL 

0115  OIST  -    S0RT(0IFFT»*2  ♦  0IFEP»*2> 

0116  IFID1ST  -  0.2)  60,61,61 

0117  61  TA(KR)  =  TT(JA) 

0118  EIKR)  =  EPIJA) 
CU9  KR  =  KR  ♦  1 
0120  IZ  =  JA 
C121  GO  TO  60 

C 

c 
c 
c 
c 

0122  62  KR  =  KR  -1 

0123  IFIKR.LE.30)  GO  TO  50 
012*  IFIKR.LE.59)   GU  TO  63 
C125  KO  =  KR/20 

Cllt  KPLCT  =  KR/KO 

CU  DO  *7   I  =  l.KELOT 

C12e  El  I)  =  El I*K0) 

0129  *7  TA I  I  )  =  TAI I*KO) 

013C  KR  =  KPLCT 

0131  GO  TO  50 

C132  63  KD  =  KR  -  30 

CliJ  KDP1  =  KD  ♦  1 

013*  DO  6*   I  =  1,KD 

0135  TAUI  =  TA(I»2) 

0136  6*  El  I  I  =  E(I»2I 

0137  DO  65   I  =  KDP1.30 
C138  TAII)  =  TAIUKC) 
C139  65  El  I  I  =  El  UKU) 
01*C  KR  =  30 

C 

C  THE  INTEGER  "NSKIF"  DETERMINES  WHICH  GRAPH  PLOTS  WILL  BE 

C  EXECUTED 

C  NSKIP  =  0:   DRAW  ONLY  THE  CREEP  CURVES 

C  NSKIP  =  1:   DRAW  ONLY  THE  CURVES  FOR  DETERMINING  C  ANE 

01*1  50  IFINSKIP.NE.il   GO  TO  818 

01*2  DO  811   I  =  1,KPTS 

Cl*3  TYSI  I  I  =  TLNZl I  I 

CI**  811  EDS! I )  =  EDZZI II 

C 

C 

C 

c 
c 

01*5  TLO  =  ABSI TYSI 1  I  ) 

01*6  [TLO  =  TLO»10. 

01*7  TLC  =  ITLO  ♦  1 

01*8  TLC  =  TLC/10. 

01*9  IFITYSI  l).(iT.O.O) 

CUC  EDLO  =  ABSIEDSt  1)  I 

C151  IEDLG  *  EDL0*1C. 

C152  EDLO  =  IEDLO  ♦  1 


THE  FOLLOWING  STATEMENTS  (DOWN  TO  STATEMENT  50)  FURTHE 
REDUCE  THE  NUMBER  OF  DATA  POINTS  TO  THIRTY  OR  LESS,  A< 
REOUIRFO  BY  SUBROUTINE  DRAW. 


THE  FOLLOWING  STATEMENTS  COMPUTE  GRAPH  SCALES  FOR  PL01 
THE  C  AND  0  DETERMINATION.   THF  NUMBER  OF  DATA  POINTS 
ARE  REDUCED  TO  THIRTY  OR  LESS. 
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0153  EDLC  »  EOLO/10. 

015*  IF(EDSIll.GT.O.O)   EDLO  '    -EOLO 

0155  FLAG    -    0.000001 

0156  02    •    ABSIITYSIKPTSI     -    TYSI1II/8.) 

0157  az   '    QZMOOOOOO. 
0156  812    az    =    0Z/10. 

0159  FLAG  *  FLAG»10. 

0160  IFtQZ  -  10.1  Hip  812,  812 

0161  813  IOZ  =  01     >     1. 

0162  ZSCAL  "  IOZ 

0163  ZSCAL  -  ZSCAL*FLAG 
016*  FLAG  =  0.000001 

0165  OZ  *  A8S( (EOS(KPTS)  -  F0S(lll/8.) 

0166  az  ■  02*1000000. 

0167  81*  8Z  =■  0Z/10. 

0168  FLAG  *  FLAG*10. 

01(5  If  (OZ  -  10.)  822.81*, 81* 

0170  822  IOZ  =  OZ  ♦  1. 

0171  YSCAL  ■  IOZ 

0172  YSCAL  -  YSCAL*FLAG 
3  73  DO  823   I  -  liKPTS 
017*  TYSII)  -  TYSlfl  ♦  TLO 
Q  75  823  EOSIII  -  EDSIII  ♦  EDLO 
6  76  MR  >  KPTS 

C  77  IFIMR.LE.30)   GO  TO  80 

0178  IFIMR.LE.59)   GO  TO  76 

C179  MO  ■=  MR/20 

0180  MPLOT  -  MR/MQ 

0181  00  75   1  -  1. MPLOT 

0182  EDSIII  -  EDSlI*MOI 

0183  75  TYSUI  -  TYS!  I*M01 
018*  MR  .  MPLOT 

0185  GO  TO  80 

0186  76  MD  =  MR  -  30 

0187  MDP1  =  MO  ♦  1 

0188  00  78   I  =  l.MD 

0189  EOS! I  )  =  E0S(I*21 

0190  78  TYSII)  =  TYS(I»2) 

0191  00  79  I  =  MDP1.30 
C192  EOS(I)  =  EDS!I»M0) 
0193  79  TYS(l)  =  TYS!I»MD) 
C19*  MR  =.  30 

0195  80  CONTINUE 

0196  LABS  ■  LABELIKCI 
C 

C  MAKE  A  PLOT  OF  THE  DATA  USED  TO  DETERMINE  C  AND  0. 
C 

0197  CALL  DRAW(MR,TYS.EDS,1,1,LA8A,KTITLE,ZSCAL.YSCAL,1,1,2,2,°,9, l.LAS 
I  TO 

C198  PRINT  815,  LASTC 

0199  815  FORMAT!//, lOX.'LASTC  =  '.121 

0200  TLCW  =  TLO 

0201  EDLOH  =  EOLO 

02C2  00  82*   1  =  l.KPTS 

0203  IUII)  =  TLNZII)  ♦  TLOW 

C2C*  82*  YUdl  =  FEDIII  ♦  EDLOW 

C 

C  SUPERIMPOSE  ON  THE  DATA  THE  FITTED  CUPVE  FOUND  FROM 

C  SUBROUTINE  LSOPOL. 

C 

0205  CALL  DRAWIKPTS,TU,YU,3,0,LABB,KTITLE,ZSCAL,YSCAL,1,1,2,2,9,9,  I,  LAS 

1TD) 

C2C6  PRINT  816,  LASTD 

0207  816  FORMAT!//, 10X, 'LASTD  =  ',  121 

0208  PRINT  817, TLO, EOLO 

0209  817  FORMAT!///, 10X, 'FOR  THE  PLOT  USEO  TO  OETERM! NF"C"  ANP  "0",  THE  SCA 

1LES  MUST  BE  CORRECTED  AS  FOLLOWS :',// .30X ,' SUBTR ACT' , F*. 2 ,' FROM  x- 
2SCALE  VALUES' ,/,30X, 'SUBTRACT', F*. 2, 'FROM  Y-SCALE  VALUFS',//) 

C 

C  IF  THE  PREVIOUS  GRAPHS  HAVE  BEEN  EXECUTED,  DO  NOT  DRAW 

C  CREEP  CURVES  BUT  GO  TO  STATEMENT  819  AND  PRINT  OUT  THE 

C  RESULTS  OF  THE  EMPIRICAL  CREEP  FIT. 

C 

021C  GO  TO  819 

0211  818  CONTINUE 
C 

C  THE  NEXT  THREE  STATEMENTS  PROVIDE  THE  PROPER  VALUE  OF  THE 

C  CALLING  ARGUMENTS  "JC  AND  "MODCUR"  SO  THE  CREEP  CURVES 

C  MAY  BE  PLOTTED. 
C 

0212  IFIJC.GT.5)   JC  =  JC  -  5 

0213  IFIL.EO.ll  MODCUR  =  1 
021*  IF!L.GT.l)  MODCUR  =  2 

0215  DO  66   I  =  l,KP 

0216  66  EHIII  ■  -  Elll 

0217  CALL  DRAW(KR,EM,TA,M0DCUR,JC,LABA, ITITLE,EPSKAL,TSKAL,C,B,2,2,8,15 


1,1, LASTA) 
PRINT 


0218  PRINT  29,  LASTA 

C219  29  FORMAT!//, 10X, 'LASTA  =  ',  121 

C22C  IFIL.EO.IXI  MODCUR  =  3 

0221  IFIL.LT.IXI  MODCUR  =  2 

0222  LABB  ■  LABELIKCI 

0223  CALL  DRAW ( 200, YGM, TG, MODCUR, 0 ,L ABB, IT  I TLE , EP SKAL , TSK AL , 0, 8 , 2 ,2, 8, 1 
15,1, LASTB) 

022*  PRINT  30,  LASTB 

0225  30  FORMAT!//, 10X, 'LASTB  =  ',  I2.////I 

C226  819  CONTINUE 

C 

C  FITTED  PARAMETERS  ARE  PRINTED  AS  OUTPUT. 

C 

0227  PRINT  27 

0228  27  F0RMATI5X, • THE  PAHAMFTERS  OF  THE  FITTED  CURVE  ARF:',//) 
C229  PRINT  33 

0230  33  F0RMATI9X, 'Al' ,12X,'A2' ,12X,'A3' ,12X, 'C , 13X,'0" ,/) 

0231  PRINT  28,  Al,A2,A3,CC.O 

0232  28  F0RMATI1H  ,  3D14.  *,  2F  1<,.*, // /  I 

0233  AA1  =  Al 
023*  AA2  =  A2 

0235  AA3  =  A3 

0236  CCC  =  CC 

0237  oeo  =  0 

C 

C  FITTED  PARAMETERS  ARE  PUNCHED  ON  CARDS  FOR  USE  IN 

C  SUBSEQUENT  PROGRAMS. 

C 

0238  PUNCH  992,  AA1 , AA2 , AA3 ,CCC , QOO 
C239  992  FORMAT(3E15.5,2F10.5) 

C2*0  DO  560   I  =  1.NPTS1 

C2*l  560  TS( II  -    TTII I 
C 

C  THE  STATEMENTS  DOWN  TO  STATEMENT  56*  ARRANGE  FOR  A 

C  SYMBOL  "P"  TO  BE  PRINTED  IN  THE  TABULATED  DATA  IF  AN 

EXPERIMENTAL  DATA  POINT  WAS  PLOTTEO  ON  THE  CREEP  CURVES. 

02*2  DO  56*   I  ■  l.NPTSl 

02*3  J  =  1 

02**  561  VAL  =  ABSITSII)  -  TAIJII 

02*5  IFIVAL  -  0.000C1I  562,562,563 

C2*6  562  PMARKII)  ■  PMKP 

02*7  GO  TO  56* 

C2*8  563  PMARKII)  =  ZILCH 

02*9  J  =  J  ♦  1 

025C  IFIJ.GT.KRI   GO  TO  56* 

0251  GO  TO  561 

0252  56*  CONTINUE 

0253  00  51   I  =  l.NPTSl 

025*  YGKII  ■  EQNlAl,  A2.A3.CC, a, TT(  I  I  I 

0255  51  ERRII)  =  IIYGllf)  -  EP! I i I *10C. I /EP! 1 1 
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C  PRINT  OUT  THE  OBSFRVEO  STRAIN,  FITTED  STRAIN,  AND 

C  PERCENT  ERROR  IN  A  POINT-BY-POINT  TABULATION. 

C 

0256  PRINT  31,  RUN 

0257  31  FORMAT! 10X, 'THE  DATA  FITTEO  AND  PLOTTED  FOR  RUN  ',A5,'  IS  LISTED  B 

1FLUW  ,///,6X. 'POINT  NO.      TIME. HOURS      ORSERVFO  STRAIN      CALC 
2  STRAIN      PCT  ERROR      «P"  IF  PLOTTED',///) 

0258  PRINT  32,  (  I  ,  TT  (  I  I  .  EP  I  I  I ,  YGH  I  I  ,  ERR  t  I  I  ,  PMARK  (  I  I  .   I  =  1.NPTS1I 

0259  32  F0RMAT(1H0,I10.F18.5,7X,F10.5,9X,F10.5,5X,F8.5,A12I 

0260  48  CONTINUE 

0261  42  CONTINUE 

0262  41  CONTINUE 
C 

C  THE  NEXT  FOUR  STATEMENTS  CAUSE  ANY  ''LEFTOVER"  DATA  SFTS  TO 

C  BE  FITTED  AND  PLOTTEO  BEFORE  ENDING  THE  PROGRAM. 

C  AT  THIS  POINT  THERE  ARE  BETWEEN  ZERO  AND  (IX  -  II  DATA 

C  SETS  NOT  YET  ANALYZED. 

C 

0263  IFINOAT  -  NDOI  53,53,52 

0264  52  IX  =  NDAT  -  NDCJ 

0265  NUC  =  NDAT 

0266  GO  TO  40 

0267  53  CONTINUE 
0266  END 


FORTRAN  IV 

G  LEVEL 

1,  MOD  1             EON 

00C1 

FUNCTION  E0N(Al,A2,A3,CC,O, Z) 
IMPLICIT  REAL'S  (A-HI.  REAL*8  IO-Z) 
1FIZ.E0.0.0I   GO  TO  704 

0CC2 

0CC3 

0004 

PUR  =  CC»(/»»0I 

0CC5 

IFIPWR  -  50. 1  702, 701,701 

CCC6 
0CC7 

701 

TERM  =  1.0 

GO  TO  703 

occe 

702 

TERM  =  1.  -  DEXP(-PWR) 

oco 

703 

EON  =  Al  ♦  A2*Z  ♦  A3*TERM 

8c°iS 

704 

GO  TO  708 
EQN  =  Al 

0C12 

7C8 

CONTINUE 

0C13 

RETURN 

0C14 

END 

OUTPUT    FROM    PROGRAM    "FITDATA" 


A2   EXTRAPOLATION,  FROM    SUBROUTINE     LSQPOL 


COEFFICIENTS  OF  Y=8 1 ♦62*X* E TC  AND  ERRORS 

Bl  11=  1.44438270-02  ERRB  =  4.065E-04  B(  21=   2. 9625 553D-05  ERRB=  5.166E-06  R( 

SIGMA=  1.0256402E-04  F  LEVEL=    0.0             SUM  SO  DEV  =    7. 5 1 384 38F-09 

CCtFFICIENTS  OF  Y= T 1 »P 1 *T2»P2*E TC  AND  ERRORS 

T(  1)=  1.67642860-02  ERRT=  3.877E-C5  T(  21=   2. 9625553D-05  ERRT=  5.166E-06  T( 

I  XIII               Fill               Yl I  I               DELYI  II            W(I) 

1  7.C000000D  CI    1.6400000D-02    1.65176160-02    1 . 1 761616E-04    1.0000000E  OC 

2  7.19667G0D  01    1  .65500000-02    1 . 557588 1D-02    2 . 5880727E-05    l.OOOOOOOE  00 

3  7.3833300D  01    1 .66000000-02    1 . 6631 180D-02    3 . 1 179792E-05    l.OOOOOOOE  OC 

4  7.7867000D  01    1 . 68500C0D-02    1. 6750680D-02   -9 . 9319601 E-05    l.OOOOOOOE  00 

5  7.8016700D  01    1.69000000-02    1 . 6755 1 150-02   -  1 . 4488466E-04    l.OOOOOOOE  00 

6  8.26C0000D  01    1 .69000000-02    1 . 68908980-02   -9. 101 8655E-06    l.OOOOOOOE  00 

7  5.40COCCCD  01    1 . 715000OD-02    1.72286290-02    7. 8629440E-05    l.OOOOOCOF  00 
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DETERMINATION    OF   PARAMETERS    "C"    A^D    "Q" 


COEFFICIENTS  OF  Y  =  B  1  »B2»X»ETC  AND  ERRORS 

bt  11=  -3.12931050-01  ERRB=  2.330E-02  B(  21=  3.98679070-01  EPRB=  8.252F-03  R( 

SIGMA'    1.4200133E-01  F  LEVEL'    0.0  SUM  SO  OEV=  1.9103099E-02 

COEFFICIENTS  OF  Y  =  1 1 »P1 *T2»P2»E TC  AND  ERRORS 

I»  11=  -1.4279344D-01  ERRT=  2.304E-02  T<  21=  3.98679070-01  ERRT*  8.252F-03  T( 

I             XI  I »              Mil  Y< I  I  OELYI II            will 


i 


-4.09414460  00   -1.80231440  00  -1.94518080  00  -1 .4286631 E-01    l.OOOOOOOE  00 

-3.68887950  00   -1.6604058D  00  -1.78361010  00  - 1. 2320423E-01    l.OOOOOOOF  00 

3  -3.40129740  00   -1.56381250  00  -1.66895710  00  -1 .0514462E-01    l.OOOOOOOF  00 

4  -3.17821380  OC   -1.47407730  00  -1.58001840  00  -1 . 0594100E-01    l.OOOOOOOF  OC 

5  -2.99573230  00   -1.44559500  00  -1.5072668D  00  -6 . 1671 790E-02    l.OOOOOOOF  OC 

6  -2.84181020  00   -1.39028630  00  -1.4459CU3D  00  -5. 56 1 49S6E-02    l.OOOOOOOE  OC 

7  -2.70800020  00   -1.36353B90  00  -1.39255400  00  -2. 901 5 1 B7E-02    l.OOOOOOOE  00 

8  -2.48494670  00   -1.28620760  00  -1.30362730  00  - 1 . 74 19655F-0?    l.OOOOOOOE  00 

9  -2.30258510  00   -1.18930850  00  -1.23092350  00  -4. 16 14994E-02    l.&OCOOOOE  00 

10  -2.14840580  00   -1.16625190  00  -1.16945550  00  -3. 20 35466E-03    l.OOOOOOOF  00 

11  -2.01492800  00   -1.12098360  00  -1.11624070  00  4. 7428943E-0 3    l.OOOOOOOF  00 

12  -1.89712000  00   -1.05545600  00  -1.06927310  00  -  1 . 3817094E-02    l.OOOOOOOE  00 

13  -1.69646750  00   -9 . 928 17 350-01  -9.89277120-01  3 .5402330 E-C3    l.OOOOCOOF  00 

14  -1.32174330  00   -8.55536100-01  -8.39882450-01  1 . 5653651 E-02    l.OOGOOCOF  00 

15  -1.04982210  00   -7. 64026 J9D-0 1  -7.31473150-01  3. 2553840E-02    l.OOOOOOOr  00 

16  -8. 36255720-01   -  7. 29 183090-C 1  -6. 4632870D-01  8 . 2854390E-02    l.OOOOOOOE  00 

17  -6.60350910-01   -6.60462630-01  -5.76199130-01  8. 4263444E-02    l.OOOOOOOE  00 

18  -2.65698820-01   -5.14083030-01  -4.18859600-01  9. 5223427E-02    l.OOOOOOOE  OC 

19  1.65325810-02   -4.06496960-01  -3.06339850-01  1 .00 1 570BE-0 1    l.OOOOOOOE  00 

20  4.16536920-01   -2.89942540-01  - 1 . 46 86650D-01  1 .4307600E-01    l.OOOOOOOF  00 

21  7.01462510-01   -1.91264420-01  -3.32726280-02  1 . 5799177F-01    l.OOOOOOOE  OC 

22  1.10416350  00    1.62578570-02  1.27275830-01  1 . 1 101 794E-01    l.OOOOOOOF  00 

23  1.39046070  00    7. 948B165D-02  2.4141652D-0 1  1 .6 1 92830F-01    l.OOOOCOOF  OC 

24  1.61277230  OC    1.86659490-01  3.30047530-01  1 .4338803E-01    l.OOOOOOOF  OC 

25  1.79453890  00    2.96263960-01  4.02514060-01  1 .06 25005F-01    l.OOCOOOOE  00 

26  2.49456470  00    5.40936280-01  6.81607640-01  1 .4067 1 31 F-01    l.OOOOOOOF  00 

27  3.15771080  OC    6. 90561770-01  9.45982150-01  2 . 5542033E-01    l.COCOOCOF  OC 

28  3.27777470  00    8. 42902470-C 1  9.93849120-01  1 . 5094662 E-01    l.OOOOOOOF  00 

29  3.39170920  00    9.18610750-01  1.03927240  00  1 . 2066162E-01    l.OOOOOOOE  00 

30  3.82173280  00    1.1166443D  00  1.21071380  00  9 .4069481 E-02    l.OOCOOOOE  00 

31  3.87258960  00    1.17866040  00  1.23098940  00  5 . 2328981 E-C2    l.OOOOOOOE  00 

32  3.9123569C  00    1.26375620  00  1.24684380  00  -  1 .69 1 2 397E-02    l.OOOOOOOE  00 

33  3.95156480  00    1.37069550  00  1.26247510  00  - 1 . 0822034E-01    l.OOCOOOOE  00 

34  3.98929330  00    1.43201240  00  1.27751670  00  -1 . 5449566E-01    l.OOOOOOOF  OC 

35  4.07073470  00    1.56382370  00  1.30998570  CO  -2. 5333794E-01    l.OOCOOCOF  00 

36  4.24849520  00    1.47937750  00  1.38085510  00  -9.852242 5E-02    l.OOOOOOOF  OC 

37  4.27620350  00    1.77565370  00  1.39190180  00  -3 . 83 751B7F-01    l.OOOOOOOE  OC 

38  4.30180980  00    1.74359600  00  1.40211050  00  -3 .4148550E-01    l.OOOOOOOE  OC 

THE  DATA  USED  TO  FIND  "C"  AND  "0"  FOR  RUN  2AP13  IS  LISTED  BELOW: 

CATA  PT.   LN  TIME     LN l-LNI 1 . -EDt I    FITTED  VALUES    PCT. ERROR 

-1.8023  -1.9452  -7.9? 

-1.6604  -1.7836  -7.42 

-1.563B  -1.6690  -6.72 

-1.4741  -1.5800  -7.19 

-1.4456  -1.5073  -4.27 

-1.3903  -1.4459  -4.00 

-1.3635  -1.3926  -2.13 

-1.2862  -1.3036  -1.35 

-1.1893  -1.2309  -3.50 

-1.1663  -1.1695  -0.27 

-1.1210  -1.1162  0.42 

-1.0555  -1.0693  -1.31 

-0.9928  -0.9893  0.36 

-0.8555  -0.8399  1.83 

-0.7640  -0.7315  4.26 

-0.7292  -0.6463  11.36 

-0.6605  -0.5762  12.76 

-0.5141  -0.4189  18.52 

-0.4065  -0.3063  24.64 

-0.2899  -0.1469  49.35 

-0.1913  -0.C333  82.60 


11 

-4.0941 

12 

-3.6889 

13 

-3.4013 

14 

-3.1782 

15 

-2.9957 

16 

-2.8418 

17 

-2.7080 

18 

-2.4849 

19 

-2.3026 

20 

-2.1484 

21 

-2.0149 

22 

-1.8971 

23 

-1.6965 

24 

-1.3217 

25 

-1.0498 

26 

-0.8363 

27 

-0.6604 

28 

-0.2657 

29 

0.0165 

30 

0.4165 

31 

0.7015 

IS 

\-.\m 

kmi 

B:UI3 

-m-.ft 

34 

1.6128 

0.1867 

C.3300 

-76.82 

35 

1.7945 

C.2963 

0.4025 

-35.86 

36 

2.4946 

0.5409 

C.6816 

-26.01 

37 

3.1577 

0.6906 

0.9460 

-36.99 

38 

3.2778 

0.8429 

0.9938 

-17.91 

39 

3.3917 

0.9186 

1.0393 

-13.14 

40 

3.8217 

1.1166 

1.2107 

-8.42 

41 

3.8726 

1.1787 

1.2310 

-4.44 

42 

3.9124 

1.2638 

1.2468 

1.34 

43 

3.9516 

1.3707 

1.2625 

7.90 

44 

3.9893 

1.4320 

1.2775 

10.79 

45 

4.0707 

1.563B 

1.3100 

16.23 

46 

4.2485 

1.4794 

1.3809 

6.66 

47 

4.2762 

1.7757 

1.3919 

21.61 

48 

4. 1018 

1.7436 

1.4021 

19.50 
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TABULATION    OF    CALCULATED     PARAMETERS    AND 
RESULTING     CURVE    FIT 


THE    PARAMETERS    OF    THE    FITTED    CURVE    AREl 

Al                                    A2                                   A3 

C 

0 

0.4957D-C2             0.29630-0*             0.9*870-02 

0.7313 

0.3987 

THE  OATA  FITTED  AND  PLOTTED  FOR  RUN  2AP13  IS  LISTED  BELOW 
POINT  NO.      TI*E, HOURS      OBSERVED  STRAIN      CALC  STRAIN       PCT  ERROR 


•P"  IF  PLOTTfD 


1 

CO 

0.00*96 

2 

0.00056 

0.00550 

3 

0.00139 

0.00560 

* 

0.00278 

0.00575 

5 

0.00*17 

0.00590 

6 

0.00556 

0.00600 

7 

0.0069* 

0.C0610 

8 

0.00833 

0.00615 

9 

0.01111 

0.00625 

10 

0.013S8 

O.C0630 

11 

0.01667 

0.006*0 

12 

0.025CO 

0.00660 

13 

0.03333 

0.00675 

14 

0.0*166 

0.00690 

15 

0. 05000 

0.00695 

0.00*96 

CO 

0.00530 

-3.611 

0.005*5 

-2.718 

0.00560 

-2.633 

0.00571 

-3.277 

0.00579 

-3.**0 

0.00587 

-3.818 

0.00593 

-3.5*2 

0.0060* 

-3.299 

0.0061* 

-2.569 

0.00622 

-2. 700 

0.006*3 

-2.6*8 

0.00659 

-2.*0  7 

0.00672 

-2.546 

0.0068* 

-1.530 

16 

0.05832 

0.00705 

17 

0.06667 

0.00710 

18 

0.08333 

0.00725 

19 

0. 10000 

0.007*5 

20 

0. 11667 

0.00750 

21 

C13333 

0.00760 

i2 

0.15000 

0.00775 

23 

0.18333 

0.00790 

2* 

0.26667 

0.00825 

25 

C.35CC0 

0.00850 

26 

0.*3333 

0.00860 

27 

0.51667 

0.00880 

28 

C76667 

0.C0925 

29 

1.01667 

0.C0960 

30 

1.5167C 

0.01000 

31 

2.01670 

0.01035 

32 

3.01670 

0.01110 

33 

4.01670 

0.01135 

3* 

5.01670 

0.01175 

35 

6.C1670 

0.01215 

36 

12.11670 

0.01310 

37 

23.51670 

0.01385 

38 

26.5167C 

0.01*30 

39 

29.7167C 

0.01*55 

*0 

45.68330 

CC1535 

*1 

48.06670 

0.01550 

*2 

50.0167C 

0.01565 

*3 

52.01670 

0.01580 

** 

54.01670 

0.C1590 

*5 

58.60000 

0.01610 

*6 

70.00000 

0.01640 

*7 

71.96670 

0.01655 

*8 

73.13330 

0.01660 

*9 

77.867C0 

0.01685 

50 

78.01670 

0.01690 

51 

82.60000 

0.01690 

52 

i*.COC(.0 

0.01715 

0.00695 

-1.423 

0.00705 

-0.760 

0.00722 

-0.475 

0.00736 

-1.173 

0.007*9 

-0.092 

0.00761 

0.14C 

0.00772 

-0.414 

0.00791 

0.109 

0.0D829 

D.502 

0.00859 

1.071 

0.00884 

2. 779 

0.0090  5 

2.RC5 

0.00=55 

3.271 

0.00993 

3.440 

0.010*9 

4.880 

0.01090 

5.290 

0.011*9 

3.479 

0.01191 

4.906 

0.01223 

4.102 

0.01250 

2.847 

0.013*9 

2.970 

0.01**2 

4.103 

0.01*59 

2.048 

0.01*76 

1.462 

0.015*7 

0.758 

0.01556 

0.3CC 

0.01563 

-0.10R 

0.01571 

-0.582 

C01578 

-0.744 

0.01595 

-0.952 

0.01634 

-0.366 

0.016*1 

-0.870 

0.016*7 

-0.794 

0.01660 

-1.476 

0.01661 

-1.739 

0.01676 

-n.855 

0.01712 

-0.172 
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APPENDIX  C 
COMPUTER  PROGRAM  FOR  PREDICTING  CREEP  BEHAVIOR  OF  COMPOSITES 
As  discussed  in  Section  7S  the  prediction  of  composite  behavior  was 
made  by  assuming  the  analytic  expression  for  creep  of  the  components  to 
be  an  equation  of  state.   Input  data  to  the  Fortran  IV  computer  program 
called  COMCUR  consists  of  the  experimental  creep  data  for  the  composite, 
the  volume  fractions  a  and  b  for  each  component,  and  the  stress  coeffi- 
cients and  exponents  for  each  component.   These  were  computed  from  stress 
correlations  of  the  creep  parameters  for  individual  components  obtained 
from  program  FITDATA.   For  a  series  of  times  ranging  over  the  testing 
period,  the  stress  in  each  component  is  computed.   This  is  done  by  a 
Newton-Raphson  iteration  since  the  unknown  cannot  be  expressed  explic- 
itly as  a  function  of  the  other  variables.   With  the  stress  distribution 
determined,  the  composite  creep  strain,  equal  to  the  strain  in  each  com- 
ponent, is  computed  using  the  equation  of  state  assumption.   Subroutine 
DRAW  is  used  to  display  the  following  output: 

(1)  Calculated  creep  curve  of  component  A  at  the  given  stress 

(2)  Calculated  creep  curve  of  component  B  at  the  given  stress 

(3)  Predicted  creep  curve  of  composite  AB  at  the  given  stress 
for  the  specified  volume  fractions  of  components  A  and  B 

(4)  Experimentally  observed  creep  data  for  composite  AB 

The  load  distribution  in  the  composite  during  the  first  200  seconds  of 
the  test  is  also  presented  in  graphical  form.   Tabulated  output  and  graphs 
for  a  typical  run  are  given  on  the  following  pages.   The  program  COMCUR 
and  the  related  function  subprograms  are  annotated  with  comment  cards. 
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COMPUTER    PROGRAM   "COMCUR' 


FORTRAN  IV  G  LEVEL  1,  HOD  1  MAIN  DATE  =  68139  08/19/41  PAGE  OOOl 

OOCl  REAL'S  ITITLEI 12 ) ,RUN,KTITLE( 121 

O0C2  REAL  KA.K8,NA,NB.MA,MB.K,M,N,LABEL( 30I.LABB  ,_.-„« 

0003  DIMENSION  fz(2C0),  STRAlJoOl,  STRB! 200  I , ECOM ( 200 ) ,  EAJ200I,  EBI20Q 
1),EMEAS!200).TMEAS(200I,FRACTA(200)  , FRACTBI 200  I , ECALC ( 200 ) , ERR ( 200 
2),EMEASM(200),EAM(200) ,EBM(200) , ECOMMI 200  I , TO < 25 ) . SA I  25 ) , SB ( 25 ) , F A 
3(25I.F8(25l,s!GG(20OI,tGI50l,EGl50l,TY(lOOI,TS(106).DA(l60l,DB(lO0 
4) 

C 

C  REAO  IN  TITLES  FOR  GRAPHICAL  OUTPUT 

C 

0004  READ  5,  (ITITLE(l),  I  ■  1.12) 
0C05  READ  5,  (KTITLE(l).  I  '  1,121 
0006  5  F0RMATI6A8I 

0CC7  REAO  6,  (LABELII).   I  *  1,3) 
00C8              6  F0RMATI3A4) 

C 

C  REAO  IN  THE  FOLLOWING  DATA: 

C  RUN  =  RUN  NUMBER  OF  EXPERIMENTAL  COMPOSITE  TEST 

C  NPTS  «  NUMBER  OF  CREEP  DATA  POINTS  FOR  THE  COMPOSITE. 

C  TMAX  =  MAX.  TIME  IN  HOURS  TO  WHICH  THE  COMPUTFD 

C  CREEP  CURVES  ARE  TO  BE  DRAWN. 

C  TSKAL  =  TIME  SCALE ( HOURS / INCH )  OF  CREEP  GRAPHS 

C  EPSKAL  =  STRAIN  SCALEIIN/IN  PER  INCH)  OF  CREEP  CURVFS 

C 

00C9  READ  7,  RUN. NPTS , TMAX ,T SKAL , E PSKAL , GL 
0010              7  F0RMAT(aB,112.4F10.0I 

C 

C  READ  IN  EXPERIMENTAL  DATA 

C  EMEAS  =  DEFORMATION  IN  MILS 

C  TMEAS  *  TIME  IN  HOURS 

OCU  REAO  8,  (EMEASdl,   I  =  l.NPTS) 

0012  8  FORMAT! 14F5.1) 

C013  READ  9,(TMEAS(II,  I  =  I, NPTS) 

0014  9  F0RMATI7F10.5) 

C 

C  REAO  IN  THE  STRESS  CORRELATED  PARAMETERS  OF  THE  EMPIRICAL 

C  CREEP  EOUATIONS  FOR  COMPONENTS  A  AND  B 

C 
0C15  REAO  10,  KA,BA,PA,OA,CA,NA,RA,MA 

0016  REAO  10,KB,B8.FB,QB,CB,NB,RB,MB 

0017  10  F0RMAT(3E15.5,5F6.4l 
C 

C  READ  IN  THE  VOLUME  FRACTIONS  (A  AND  B)  OF  EACH  COMPONENT 

C  AND  THE  STRESS  LEVEL  OF  THE  EXPERIMENTAL  TEST  (  STP  ) 

C 
0C18  READ  11,  A.B 

CC1S  11  FORMAT(2Fl6.0) 

0020  REAO  12,  STR 

0021  12  FORMAT(FIO.O) 
C 

C  GENERATE  200  TIME  INTERVALS  AT  WHICH  TO  PREDICT  THF  CREFP 

C  CURVE. 

C 

OC22  TINC  -    TMAX/200. 

0023  TZ(1)  =  TINC 

0C24  00  31  I    -    2.2C0 

0025  31  TZ( I )  =  TZ( I-ll  ♦  TINC 

0026  J  =  1 
C 

C  ESTIMATE  THE  STRESS  IN  COMPONENT  B  AT  THE  START  OF  THE 

C  TEST 

C 

0027  GUESS  =  STR»B 
0C28  00  32   I  =  1.2C0 

C 

C  IN  THE  FUNCTION  SUBPROGRAM  "STRESS",  COMPUTE  THE  STRESS 

C  IN  COMPONENT  B  AT  THE  200  TIME  INTERVALS  (STRB).   USE  THE 

C  LAST  COMPUTED  VALUE  OF  STRB  AS  THE  NEW  GUESS  AT  THE  NEXT 

C  TIME  INTERVAL. 

C 

0C29  STRBII)  «  STRESS(GUFSS,TZ( I ) , KA , KB, PA , PB , 8A , BB.NA ,NB, RA, RB, MA ,M8, 
1CA,CB,0A,0B,A,B,STRI 

0030  IFISTRBI I). EO. 50000. )   GO  TO  774 

0031  J  *  J  ♦  1 

OC32  32  GUESS  =  STPBI I ) 

CC33  DO  33   I  =  1,200 

C 

C  CALCULATE  THE  STRESS  IN  COMPONENT  A  (STRAI. 

C 
0C34  STRAUI  =  (STR  -  B*STRB(II)/A 

C 

C  USE  FUNCTION  SUBPROGRAM  "STRAIN"  TO  COMPUTE  THE  FOLLOWING: 

C  ECOM  =  PREDICTED  CREEP  STRAIN  IN  COMPOSITE  WHEN 

C  LOADED  AT  STRESS  LEVEL  EQUAL  "STR" 

C  EA  =  COMPUTED  STRAIN  IN  COMPONENT  A  ALONE  WHEN 

C  LOADED  AT  STRESS  EQUAL  TO  "STR". 

C  EB  =  COMPUTED  STRAIN  IN  COMPONENT  B  ALONE  WHEN  LOADED 

C  AT  STRESS  EQUAL  TO  "STR". 

C 
0C35  ECOM(I)  =  STRAINISTRB! I ) ,TZ( I ) , KB, PB , BB , N8 ,RB , MB ,CB , OB  I 

CC36  EA(l)  =  STRAINISTR.TZI I ) , K A , PA, BA , NA , RA , MA ,C A , Q4 I 

CG37  EB(I)  =  STRAINISTR.TZ I  I ) , KB, PB, BB, NB , RB , MB.CB, OB  I 

C 

C  CHANGE  THE  SIGNS  FOR  CLARITY  ON  GRAPHICAL  PRESENTATION. 

C 
CC36  EAMI  I  I  =    -EA!  I) 

0C39  EBM( I )  =  -EB( I ) 

CC4C  33  ECOMMI I  I  =  -ECCM(I | 

C 

C  CHANGE  EXPERIMENTAL  DEFORMATION  IN  MRS  TO  STRAIN  (IN/IN). 

C 
CC41  DO  34   I  =  l.NPTS 

0042  EMEASdl  =  EMEASII  1/(1000. *GLI 

0C43  34  EMEASMdl  =  -EMEASdl 

C 

C  THE  NEXT  TWELVE  STATEMENTS  USE  SUBROUTINE  DRAW  TO  PLOT  THE 

C  COMPUTED  CREEP  CURVE  FOR  THE  INDIVIDUAL  COMPONENTS  AND 

C  FOR  THE  COMPOSITE. 

C 
CC44  LABB  =  LABEL(l) 

0045  CALL  DR AW (200.E AM,TZ,1,0,L ABB, I T I TL F , EPSKAL , T SKAL, 0,8, 2,2,8, 15, 1, 

1LASTA) 
0C46  PRINT  40.LASTA 

0C47  40  FORMAT! //,10X. ■  LASTA  =  ',12,//) 

0048  LABB  =  LABEL<2) 

0049  CALL  DRAW! 200 , FBM, T Z , 2 ,0 , LABB, I T ITLE .EPSKAL, TSKALiO, 8,2 , 2, 8 , 15, 1, 
1LASTB) 

CC5C  PRINT  41.LASTB 

0051  41  FORMAT! //,10X,  •  LASTB  =  ',12,//) 

CC52  LABB  =  LABELI3) 

00  5  3  CALL  DRAW! 200, ECOMM.TZ, 2,0, LABB, I T I TLE , EPSKAL , TSKAL , 0 , 8 , 2 , 2 , 8 , I  5, 1 

l.LASTM) 
CC54  PRINT  42.LASTM 

0C55  42  FORMAT!//, 10X, •  LASTM  =  ',12, //I 
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THE  NEXT  TWENTY  STATEMENTS  ARE  USED  TO  SELECT  30  EXPERI- 
MENTAL OATA  POINTS,  THE  MAXIMUM  ALLOWED  IN  SUBROUTINE 


ORAW.   THE  SELECTED  EXPERIMENTAL  POINTS  ARE  SUPERIMPOSED 
--OICT-- 


C  ON  THE  PREDICTED  CREEP  CURVE. 

C 

"FINPTS.GT.59I   GO  TO  T76 


fF(NPTS.GT,30l   GO  TO  50 
IPTS  »  NPTS 

■  IPTS 


0059  00  4?   I  «  l.IPTS 
00*0  EGIII  ■  ENEASHWI 

0061  49  TGI ii  •  TMEASI I  I 

0062  GO  TO  53 

0063  50  JF  »  NPTS  -  30 
0C64  00  51   I  «  l.JF 

0065  EGIII    '    EMEASMI2*II 

0066  51     TGUI    »    THEASI2*!) 

0067  JF2   -    I2*JF >    ♦    1 

O^t*  EGlf-JF      •   efwIlSHJ  1 1 

0070  52    TG(|-JFj    -    TMEASIt) 

55 7 5  CALL    0RAW(IPTS,EG,TG,3, J,RUN,ITITLE,EPSXAL,TSKAL,0,8,2,2,8,15,1,LA 

0074  PRINT  43.1AST0 

0075  43  FORMAT!//, 10X,'  LASTO  -  ',12, //I 
C 

C  GENERATE  100  TIMES  1TYI  COVERING  THE  FIRST  200  SECONDS 

C  OF  THE  TEST.   USE  FUNCTION  SUBPROGRAM  "STRESS"  TO  FIND 

C  THE  CORRESPONDING  STUESSES  1DBI  IN  COMPONENT  B. 
C 

0076  FINC  -  2./3600. 

0077  TYI11  =  FINC 
OC78  00  70   I  *  2,100 

0079  70  TY(I)  =  TY(I-l)  ♦  FINC 

0060  J  -  1 

0081  GUESS  ■  STR»B 

0C82  00  73   I  *    1,100 

0083  DBIII  ■  STRESSIGUESS.TYI I  I , KA.KB , PA , PB, BA, BB, NA, NR, RA, RB, MA , MB, CA, 
lCB.OA.OB.A.B.STR) 

oce<i  ificbIii.eo.soooo.)     go  to  774 

0085  J   *    J    ♦    I 

0086  73  GUESS  *  OBI  I) 
C 

C  CALCULATE  THE  STRESS  IN  COMPONENT  A  (DAI  AND  THF 

C  LOAD  FRACTIONS: 

C  FRACTA  =  FRACTION  OF  THE  LOAD  CARRIED  RY  COMPONENT  A. 

C  FRACTB  =  FRACTION  OF  THE  LOAD  CARRIED  BY  COMPONENT  B. 

C 

0087  00  74   I  «  1,100 

0C88  DAIII  »  ISTR  -  B*DB(III/A 

0089  FRACTAdl  =  -IDA  I  I  l*A  l/STR 

0090  FRACTBII)  =  -i DBI I )»B ) /STR 
0C91  7*  ISIII  *    TY(II»3600. 

0092  LABB  =  LABELI1) 

0093  KTITLEI12)  *  RUN 
C 

C  USE  SUBROUTINE  ORAW  TO  PLOT  THE  LOAD  DISTRIBUTION  DURING 

C  THE  FIRST  200  SECONDS. 

C 

0094  CALL  0RAWI100, FRACTA, TS, 1,0, LABB, KTITLE, 0.2, 20.  ,0, 6  ,2 , 2 , 6, 1 5  , 1  ,L A 
1STA) 

0095  PRINT  40.LASTA 

0096  LABB  =  LABELI2I 

0C97  CALL  DRAWI100, FRACTB, TS. 3,0, LABB, KTITLE, 0.2, 20.  , 0 , 6 ,2 , 2,6, 15 , 1 ,LA 

1STBI 

CCS8  PRINT  41.LASTB 

0C99  DO  75   I  =  1,100 

0100  FRACTAdl  -  -FPACTAII) 

0101  75  FRACTBII)  =  -FRACTBII! 

0102  J  *  1 
C 

C  CALCULATE  THE  STRESS  ISIGG)  IN  COMPONENT  B  (USING  FUNCTION 

C  SUBPROGRAM  "STRESS")  CORRFSPONOI NG  TO  THE  EXPERIMENTAL 

C  TIMES  (TMEASI.   USE  FUNCTION  SUBPROGRAM  "STRAIN"  TO  FIND 

C  THE  PREDICTED  STRAIN  IECALC)  CORRESPONDING  TO  THE  DATA 

C  POINTS  SO  A  TABULATEO  COMPARISON  MAY  BE  MADE. 
C 

0103  GUESS  =  STRB(l) 

0104  00  36   I  •  l.NPTS 

01C5  SIGGII)  =  STRESSIGUESS, TMEASI I ) ,KA, KB , PA , PR, BA, BB , NA , NB , RA , RB ,MA , 

IMB.CA.CB.OA.OB.A.B.STR) 

01C6  IFISIGGI I). EC. 50000.)   GO  TO  774 

01C7  ECALCII)  =  STRAINISIGGI I ), TMEASI I > , K8 . PB , BB , NH ,RS, MB.CB.OB I 

01C8  J  =  J*l 

0109  36  GUESS  =  SIGGI I  I 
C 

C  COMPUTE  THE  PERCENT  ERROR  IN  PREOICTED  STRAIN 
C 

0110  00  37   I  =  l.NPTS 

CI  11  37  ERRIII  =  (1EMEASII)  -  ECALC I  I) I /EME AS  I  I ) )«100. 
C 

C  FIND  THE  LOAD  FRACTIONS  (FA, FBI  ANO  THE  STRESSES  AT  25 

C  SELECTED  TIMES  (TO)  COVERING  THE  DURATION  OF  THE  TEST. 

0112  00  38   I  ■  1,25 

0113  FA(I)  =  (STRA1I»8-7I*A)/STR 

0114  FBII)  =  ISTRBI l*8-7)»8l/STR 

0115  TO(I)  =  T2(I»8-7) 

0116  SAID  ■  STRAII*8-7) 
C117  38  SBI I)  =  STRB(I*R-7) 

C 

C  PRINT  OUT  THE  COMPUTED  RESULTS  IN  TABULAR  FORM. 

C 

0118  PRINT  20,  RUN 

0119  20  FORMAT!/////, 10X,'  COMPARISON  OF  THE  CALCULATED  AND  EXPERIMENTAL  V 

1ALUES  OF  CREEP  STRA I N • , / , 10X, •  FOR  COMPOSITE  RUN',A8,'  IS  GIVEN  RE 
2L0W: • ,///! 

0120  PRINT  21 

0U1  21  F0RMATI10X,'  TIME, HOURS        EXPTL.  STRAIN          CALC.  STRAIN 

1  PCT.  ERROR'. //I 

0122  PRINT  22.ITMEAS! I ) , EME AS  I  I  I , ECALC I  I l.ERRII ),  I  =  l.NPTS) 

0123  ZZ    FORMAT! 1H0.F20. 5 ,F 18. 6, F20.6.F15.2 ) 

0124  PRINT  23 

0125  23  FORMAT!/////, 10X,'  THE  CALCULATED  STRESSES  CARRIED  BY  EACH  MATERIA 

1L  OF  THE  COMPOSITE  ANO  THE  CORR ESPOND ING • , / , 10X, •  LOAO  FRACTIONS  A 
2RE  GIVEN  BELOW  FOR  25  SELECTED  TIMES  COVERING  THE  DURATION  OF  THE 
3TEST',///! 

0126  PRINT  24 

0127  24  F0RMATI7X,'  TIME,  STRESS  IN          STRESS  IN       LflAO  FRACT 

HON     LOAO  FRACTION' ,/,7X,  '  HOURS        MATERIAL  "A"      MATERIAL 
2"B"     CARRIED  BY  "A"    CARRIED  BY  "B"',//) 

0128  PRINT  25, ITOIII.SA! II.SB! I I.FAI I I.FBII),   I  =  1,25) 

0129  25  F0RMAT(1H0,F13.5,F15.I,F17.I,F16.5,F17.3) 
C13C  PRINT  26 

OUl  26  FORMAT!/////,  10X,'  THE  LOAO  FRACTIONS  AND  STRESSES  CARRIED  BY  EACH 
l',/,10X,'  COMPONENT  DURING  THE  INITIAL  PART  OF  THE  TEST  ARE  LISTED 

2  BELOW:',///) 
01J2  PRINT  27 

0133  27  F0RMATI7X,'  TIME,  STRESS  IN          STRESS  IN       LOAD  FRACT 

HON     LOAD  FRACTION'  ,/,7X.  '  SECONDS     MATERIAL  "A"      MATERIAL 
2"B"     CARRIED  BY  "A"    CARRIEO  BY  "B"'.//l 

0134  PRINT  28, (TSII I.DAI  I), DBI I), FRACTA!  I I, FRACTB! 1 ),   I  =  1,1001 

0135  28  FORMAT! 1H0 ,F 1 3.0 , F 15. 1 , Fl 7. 1 ,F 1 6. 3, F 17. 3 ) 

0136  GO  TO  776 

0137  774  PRINT  775, J 

C138  775  F0RMATI////.20X,'  CONVERGENCE  FAILURE,  J  =',I3,///I 

0139  776  CONTINUE 

C14C  END 
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FORTRAN 
GOOl 


IV  G  LEVEL  1,  MOO  1 


0CC2 
O0C3 
0CC4 
0005 
0CC6 

00C7 

ooce 

0CC9 
0010 

8811 

CC13 
0C14 
0C15 
0C16 
0C17 
0C18 
0C19 


G02C 
OCtl 
0C22 


CC23 
0024 
0025 
0026 
C027 
0C28 
0C29 


FUNCTION  STRESS(SIGB,T,KA,KB,PA,PB,BA,BB,NA,NB.,RA,RB,MA,MB,CA,CB, 

THIS  SUBPROGRAM  USES  A  NEWTON-R APHSON  ITFRATION  TO 
COMPUTE  THE  STRESS  IN  COMPONENT  B  (SIGB)  FOR  A  GIVEN 
TIME  <TI.   IN  THE  FUNCTION  (FUNCTI  THE  ONLY  UNKNOWN  IS 
SIGB.   THE  FUNCTION  IS  FORMED  BY  EQUATING  THE  STRAIN  IN 
EACH  COMPONENT.   THE  STRAIN  IN  EACH  COMPONENT  IS  GIVEN  BY 
THE  PREVIOUSLY  DETERMINED  EMPIRICAL  CREFP  EOUATION, 
WHOSF  PARAMETERS  ARE  READ  IN  AS  ARGUMENTS  OF  THE 
SUBPROGRAM.   THE  FUNCTION  IS  FORMED  BY: 

FUNCT  =  STRAIN  IN  A  -  STRAIN  IN  B  =  0. 
THE  DERIVATIVE  (DERIVI  OF  THE  FUNCTION  IS  TAKEN 
ANALYTICALLY. 

1QA,QB,A,B,STR) 
REAL  KA,K8,NA,NB,MA,M8 
L  -    0 
101  SIGBOL  ■  SIGB 

SIGA  =  (STR  -  6*SIGB)/A 

TERM1    =    KA»(SIGA*»NA)    -    KB«  (  S  IGB»*NB  I     ♦    BA*(  S  I  GA**R  A  )  *T    -    BBHSIGB 
1*»RB)*T 
PWRA    =    CA»(T«»OA) 
IFIPWRA-50.1    82.82,81 

81  GGA    •    1.0 
GO    TO    83 

82  GGA    -    1,    -    EXP(-PWRA) 

83  PWRB    *    CB»(T»«OB) 
IFIPMRB-50.1     85,85,84 

8*    GGB    =    1.0 
GO    TO    86 

85  GGB    =    1.    -    EXP(-PWRB) 

86  FUNCT  =  TERM1  ♦  PA» ( S I GA*»MA ) »GGA  -  PB* I S IGR**MB I *GGB 
H  =  -B/A 

OERIV  =  H»(NA«KA*ISIGA«»INA-1.I I  ♦  RA*BA* ( S I GA«* < RA- 1 . I ) *T  *  MA*PA 
l»(SIGA»*IMA-l. ) )*GGA)  -  (NB*KB*( SIGB*»(NB-1. ) )  ♦  RR*BB* I S IGB**( RB- 
21.))*T  ♦  MB»PB*( SIGB»»(MB-1. I I*GGB1 

SIGB  =  SIGB  -  (FUNCT/DERIV) 

DIFF  =  ABSISIGB  -  SIGBOLI 

L  -  L»l 

IF  THE  ROOT  OF  FUNCT  (=  SIGB)  IS  NOT  FOUNO  IN  TWENTY 
ITERATIONS,  ASSUME  CONVERGENCE  FAILURE. 

IFIL.GT.20t  GO  TO  94 
IFIOIFF  -  0.5)  93,93,101 

93  STRESS  =  SIGB 
RETURN 

94  STRESS  =  50000. 
RETURN 

ENC 


FGRTRAN  IV  G  LEVEL  1,  MOD  1 


DATE  =  68139 


CCC3 
0004 
O0C5 
C0C6 
CCC7 
0C08 
CCC9 
CC10 


FUNCTION  STRAIMS,2,K,P,B,N,R,M,C0) 
REAL  K,P,N 

THIS  SUBPROGRAM  USES  THE  CREEP  EOUATION,  WHOSE  PARAMETEPS 
ARE  READ  IN  AS  SUBPROGRAM  ARGUMENTS,  TO  COMPUTE  THE  STRAIN 
FOR  A  GIVEN  STRESS(S)  AND  A  GIVEN  TlMEIZ). 

POWER  =  C*IZ**0( 

IFIPOWER  -50.)  62,62,61 

PROD  =  1.0 

GO  TO  63 

PROD  =    l.-EXP(-POWER) 

STRAIN  =  K«(S**N)  ♦  B*(S»*R)*Z  ♦  P» < S«*M I IPROD 

RETURN 

END 
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TABULATED   OUTPUT  FROM 
PROGRAM   COMCUR 
FOR    PE/PP    COMPOSITE 
STRESS  *  IOOOPSI 


COMPARISON  OF  THE  CALCULATED  AND  EXPERIMENTAL  VALUES  OF  CREEP  STRAIN 
FOR  COMPOSITE  RUN     AA11  IS  GIVEN  BELOW: 


TIME. HOURS 

0.00056 


EXPTL.  STRAIN 


CALC.  STRAIN 


PCT.  ERROR 


0.00139 

0.006200 

0.005013 

19.15 

0.00278 

0.006*00 

0.00527* 

17.60 

C. 00*17 

0.006650 

0.005*56 

17.96 

0.00556 

0.006800 

0.005599 

17.66 

0.00694 

0.006900 

0.005719 

17.12 

0.00833 

0.007000 

0.005823 

16.81 

0.01111 

0.007150 

0.006001 

16.08 

0.01388 

0.007300 

0.0061*8 

15.78 

0.01667 

0.007*00 

0.006277 

15.18 

0.02222 

0.007550 

0.006*92 

1*.01 

0.02777 

0.007650 

0.006672 

12.79 

0.03330 

0.007750 

0.006826 

11.92 

O.OSOOO 

0.007900 

0.007199 

8.R7 

0.08333 

0.008350 

0.00772* 

7.50 

0.13333 

0. 008650 

0.00826* 

*.*7 

0.21667 

0.009050 

0.008879 

1.89 

0.30000 

0.009350 

0.009323 

0.29 

0.38333 

0.010050 

0.00967* 

3.7* 

0.63333 

0.010500 

0.010*30 

0.66 

1.80000 

0.011800 

0.012122 

-2.73 

3.30000 

0.012900 

0.013128 

-1.77 

3.80000 

0.013050 

0.013361 

-2.38 

*. 80000 

0.013500 

0.0137*1 

-1.79 

6.03330 

0.013900 

0.01*108 

-l.*9 

10.08300 

0.01*600 

0.01*90* 

-2.09 

21.6*999 

0.015900 

0.016035 

-0.85 

23.6*999 

0.015900 

0.016165 

-1.67 

27.16699 

0.016250 

0.016371 

-0.7* 

28.79999 

0.016*00 

0.016*59 

-0.36 

30.0*999 

0.016*00 

0.016523 

-0.75 

3*. 29999 

0.016650 

0.016727 

-0.*6 

*5. 63300 

0.017200 

0.017193 

0.0* 

*8. 01669 

0.017300 

0.017281 

0.11 

50.01669 

0.017*00 

0.017353 

0.27 

51.86699 

0.017500 

0.017*18 

0.*7 

53.8*999 

0.017550 

0.017*87 

0.36 

58.20000 

0.017600 

0.017633 

-0.19 

69.93300 

0.018000 

0.018003 

-0.02 

71.83299 

0.018000 

0.018061 

-0.3* 

73.59999 

0.018050 

0.01811* 

-0.36 

75.93300 

0.018100 

0.01818* 

-0.*6 

77.96700 

0.018200 

0.0182*3 

-0.2* 

61.66699 

0.018250 

0.018351 

-0.55 

93.*3300 

0.018550 

0.018685 

-0.73 
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THE  CALCULATED  STRESSES  CARRIED  BY  EACH  MATERIAL  OF  THE  COMPOSITE  AND  THE  CORRF  SPDND I  Mr, 
LOAD  FRACTIONS  ARE  GIVEN  BELOW  FOR  25  SELECTED  TIMES  COVERING  THE  DURATION  OF  THF  TEST 


TIME, 

HOUHS 

STRESS  IN 
MATERIAL  "A" 

STRESS  IN 
MATERIAL  "B" 

LOAD  FRACTION 
CARRIED  BY  MA" 

LOAD  FRACTION 
CARRIED  RY  "fl" 

c.soooc 

852.6 

1141.0 

0.417 

0.583 

4. 50000 

824.9 

1167.5 

0.40  3 

0.597 

6.5OG00 

821.2 

1171.1 

0.402 

C.598 

12.50000 

819.9 

1172.4 

0.401 

0.599 

16.50000 

819.3 

1172.9 

0.401 

C.599 

2C. 50000 

819.1 

1173.1 

0.401 

0.599 

24.50000 

818.9 

1173.3 

0.400 

0.600 

28.50000 

818.9 

1173.3 

0.400 

0.600 

32.50000 

818.8 

1173.4 

0.400 

0.600 

36.50000 

818.  T 

1173.5 

0.400 

0.600 

A  = 

PE 

40.50000 

818.6 

1173.6 

0.400 

0.600 

44.50000 

818.5 

1173.7 

0.400 

C.600 

B  = 

PP 

48.50000 

818.4 

1173.8 

0.400 

C.600 

52.50000 

818.3 

1173.9 

0.400 

0.600 

56.50000 

818.1 

1174.0 

0.400 

0.600 

60.50000 

818.0 

1174.2 

0.400 

0.600 

64.50000 

817.8 

1174.3 

0.400 

0.600 

68.50000 

817.6 

1174.5 

0.400 

0.600 

72.5000C 

817.4 

1174.7 

0.400 

0.600 

76.50000 

817.2 

1174.9 

0.400 

C.600 

80.5000C 

817.0 

1175.1 

0.400 

0.600 

S4.5000C 

816.8 

1175.3 

0.399 

0.601 

88.5000C 

816.6 

1175.5 

0.399 

0.601 

92.50000 

816.4 

1175.7 

0.399 

0.601 

56. 50000 

816.1 

1176.0 

0.399 

0.601 

THE  LOAD  FRACTIONS  AND  STRESSES  CARRIED  BY  EACH 

COMPONENT  DURING  THE  INITIAL  PART  OF  THE  TEST  ARE  LISTED  BELOW: 


TIME, 
SECONOS 

STRESS  IN 
MATERIAL  "A" 

STRESS  IN 
MATERIAL  »B" 

LOAD  FRACTION 
CARRIED  BY  "A" 

LOAD  FRACTION 
CARRIED  BY  "fl 

2. 

939.9 

1057.5 

0.460 

0.540 

4. 

936.5 

1060.7 

0.458 

0.542 

6. 

933.9 

1063.2 

0.457 

0.543 

a. 

931.8 

1065.2 

0.456 

0.544 

10. 

930.0 

1067.0 

0.455 

0.545 

12. 

928.4 

1068.5 

0.454 

0.546 

14. 

927.0 

1C69.9 

0.453 

0.547 

16. 

925.7 

1071.1 

0.453 

0.547 

18. 

924.5 

1072.3 

0.452 

0.548 

20. 

923.3 

1073.4 

0.452 

0.548 

22. 

922.3 

1074.3 

0.451 

0.549 

24. 

921.3 

1075.3 

0.451 

0.549 

26. 

920.4 

1076.1 

0.450 

0.550 

28. 

919.6 

1077.0 

0.450 

0.550 

30. 

SIB. 7 

1C77.8 

0.449 

0.551 

32. 

918.0 

1C78.5 

0.449 

0.551 

34. 

917.2 

1079.2 

0.449 

0.551 

36. 

916.5 

1079.9 

0.448 

0.552 

38. 

915.8 

1C80.6 

0.448 

0.552 

40. 

915.2 

1C81.2 

0.448 

0.552 

42. 

914.5 

1081.8 

0.447 

0.553 

44. 

913.9 

1082.4 

0.447 

0.553 

46. 

913.3 

1083.0 

0.447 

0.553 

48. 

912.7 

1C83.5 

0.446 

0.554 

50. 

912.2 

1084.0 

0.446 

0.554 

52. 

911.7 

1084.5 

0.446 

0.554 

54. 

911.1 

1C85.0 

0.446 

0.554 

56. 

910.6 

1085.5 

0.445 

0.555 

58. 

910.1 

1CR6.0 

0.445 

0.555 

60. 

909.7 

1086.4 

0.445 

0.555 
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62. 

909.2 

6*. 

908.8 

66. 

908.3 

68. 

907.9 

70. 

907.5 

72. 

907.0 

7*. 

906.6 

76. 

906.3 

78. 

905.9 

80. 

905.5 

82. 

905.1 

84. 

904.8 

86. 

904.4 

88. 

904.1 

90. 

903.7 

92. 

903.4 

94. 

903.1 

96. 

902.7 

98. 

902.4 

100. 

902.1 

102. 

901.8 

104. 

901.5 

106. 

901.2 

108. 

900.9 

110. 

900.6 

112. 

900.3 

114. 

900.1 

116. 

899.8 

118. 

899.5 

120. 

899.2 

122. 

899.0 

124. 

898.7 

126. 

898.5 

128. 

898.2 

130. 

898.0 

132. 

897.7 

134. 

897.5 

136. 

897.2 

138. 

897.0 

140. 

896.8 

142. 

896.5 

144. 

896.3 

146. 

896.1 

1*8. 

895.9 

150. 

895.6 

152. 

895.4 

1S4. 

895.2 

1086.9 

0.445 

1087.3 

0.444 

1087.7 

0.444 

1088.2 

0.444 

1088.6 

0.444 

1088.9 

0.444 

1089.} 

0.441 

1089.7 

0.443 

1090.1 

0.443 

1090.4 

0.443 

1090.8 

0.443 

1091.1 

0.442 

1091.5 

0.442 

1091.8 

0.442 

1092.1 

0.442 

1092.5 

0.442 

1092.8 

0.442 

1093.1 

0.441 

1093.4 

0.441 

1093.7 

0.441 

1094.0 

0.441 

1094.3 

0.441 

1094.5 

0.441 

1094.8 

0.441 

1095.1 

0.440 

1095.4 

0.440 

1095.6 

0.440 

1095.9 

0.440 

1096.2 

0.440 

1096.4 

0.440 

1096.7 

0.440 

1096.9 

0.439 

1097.2 

0.439 

1097.4 

0.439 

1C97.6 

0.439 

1C97.9 

0.439 

1C98.1 

0.439 

1098.3 

0.439 

1098.6 

0.439 

1098.8 

0.439 

1099.0 

0.438 

1099.2 

0.438 

1099.4 

0.438 

1C99.7 

0.438 

1099.9 

0.438 

1100.1 

0.438 

1100.3 

0.438 

156. 

895. 0 

1100.5 

0.438 

158. 

894.8 

1100.7 

0.438 

loO. 

894.6 

1100.9 

0.437 

162. 

894.4 

1101.1 

0.437 

164. 

894.2 

1101.3 

0.437 

16b. 

894.0 

1101.5 

0.437 

168. 

893.8 

1101.6 

0.437 

170. 

893.6 

1101.8 

0.437 

172. 

893.4 

1102.0 

0.437 

174. 

893.2 

1102.2 

0.437 

176. 

893.0 

1102.4 

0.437 

178. 

892.8 

1102.6 

0.437 

180. 

892.6 

1102.7 

0.436 

182. 

892.4 

1102.9 

0.436 

184. 

892.3 

1103.1 

0.436 

186. 

892.1 

1103.3 

0.436 

188. 

891.9 

1103.4 

0.436 

190. 

891.7 

1103.6 

0.436 

192. 

891.5 

1103.8 

0.436 

194. 

891.4 

1104.0 

0.436 

196. 

891.2 

1104.1 

0.436 

H8. 

891.0 

1104.3 

0.436 

200. 

890.9 

1104.4 

0.436 

0.555 
0.556 
0.*56 
0.556 
0.556 
0.556 
0.557 
0.557 
0.557 
0.557 
0.557 
0.558 
0.558 
0.558 
0.558 
0.558 
0.558 
0.559 
0.559 
0.559 
0.559 
0.559 
0.559 
0.559 
0.560 
0.560 
0.560 
0.560 
0.560 
0.560 
0.560 
0.561 
0.561 
0.561 
0.561 
0.561 
0.561 
0.561 
0.561 
0.561 
0.562 
0.562 
0.562 
0.562 
0.562 
0.562 
0.562 


0.562 
0.562 
0.563 
0.563 
0.563 
0.563 
0.563 
0.563 
0.563 
0.563 
0.563 
0.563 
0.564 
0.564 
C.564 
0.564 
0.564 
0.564 
0.564 
0.564 
0.564 
0.564 
0.564 
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APPENDIX  D 
PREDICTION  OF  STRESS  RELAXATION  BEHAVIOR  FROM  EMPIRICAL  CREEP  EQUATION 

A  stress  relaxation  test  is  performed  by  rapidly  deforming  a  material 
to  a  predetermined  strain  then  observing  the  decrease  in  stress  as  a  func- 
tion of  time  while  this  strain  is  maintained.   A  series  of  such  tests  were 
run  on  the  Instron  testing  machine  for  the  plastics  used  in  this  study. 
The  empirical  equation  developed  for  creep  behavior  was  used  to  predict 
the  relaxation  curve. 

-ctQ 


(D.l) 


£=KdN+  &crRt   +  ?<S*(\-£^  ) 


Again  the  equation  of  state  assumption  was  made  for  predicting  relaxation 
behavior.   Thus,  for  a  given  strain,  and  a  series  of  times,  the  corres- 
ponding stress  may  be  found  by  implicitly  solving  equation  (D.l).   This 
was  done  by  a  Fortran  IV  computer  program  called  program  RELAX.   The  input 
data  consisted  of  the  creep  parameters  previously  determined  and  the  ini- 
tial strain.   A  Newton-Raphson  iteration  was  used  to  solve  equation  (D.l) 
for  stress  as  a  function  of  time,  and  these  values  were  tabulated  as  out- 
put.  A  plotting  subroutine  available  at  the  NPGS  Computer  Facility  (SUB- 
ROUTINE DRAW)  was  used  to  display  the  predicted  relaxation  curves.   Fig- 
ure 32  illustrates  a  computed  relaxation  curve  compared  with  the  experi- 
mental relaxation  curve  during  the  initial  part  of  the  test.   Figure  33 
shows  the  predicted  curve  and  experimental  data  over  the  entire  time  frame 
of  the  experiment.   These  results  are  typical  of  all  runs;  the  same  general 
trend  was  observed  for  each  plastic,  even  at  low  stress  levels.   The 
empirical  creep  equation  predicts  a  greater  decrease  in  stress  than  that 
observed  experimentally.   Apparently  there  is  a  structure  factor  influence 
creating  irreversibilities  that  are  unaccounted  for  in  the  creep  equation. 
The  equation  of  state  assumption  is  not  valid  --  the  stress  level  depends 
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on  the  strain  history. 

An  interesting  observation  can  be  made  from  the  experimental  data 

in  Figure  32.   During  the  first  ten  seconds  of  the  test,  the  stress  fell 
off  rapidly  in  exponential  fashion.   This  indicates  the  presence  of  mole- 
cular mechanisms  with  relaxation  times  on  the  order  of  several  seconds. 
The  corresponding  retardation  time  in  the  creep  equation  is  approximately 
one  hour;  the  predicted  curve  does  not  show  this  sudden  decrease  at  the 
start  of  the  test.   It  can  be  concluded,  therefore,  that  the  "instantaneous 
strain",  parameter  A1  in  the  creep  equation,  includes  a  portion  of  the  de- 
layed strain  with  very  short  time  constants.   The  time  scale  of  the  creep 
test  and  the  method  of  extrapolation  to  determine  parameter  A-  prevents 
the  resolution  of  this  rapid  but  delayed  elastic  response.   From  Appendix 
A,  Table  3,  it  is  seen  that  the  stress  exponent  N  for  the  instantaneous 
strain  is  slightly  greater  than  1.   If  the  delayed  strain  component  com- 
bined in  the  term  \\C5       could  be  removed,  this  exponent  would  decrease. 
Although  not  a  necessary  consequence  of  the  theory  of  rubber-like  elasti- 
city [2],  it  appears  that  the  exponent  may  approach  a  value  of  1,  so  that 
under  instantaneous  loading  the  plastics  behave  as  ideal  Hookean  solids. 

A  modified  equation  was  used  with  program  RELAX,  separating  the  para- 
meter A1  into  an  instantaneous  response  and  a  rapid  delayed  elastic  re- 
sponse: 

(d.2)  e  ■=--    K,(T      +   K-z  G   (I  -  X    ) 

where 

(d.3)  K  *=  K.  +   K?. 

The  stress  coefficients  K.  and  K„  were  estimated  from  the  experiment- 

1      I 

al  relaxation  data  and  the  time  constant  '7a.  taken  to  be  10  seconds. 

The  predicted  relaxation  curve  then  matched  the  data  more  closely  at  short 

times,  but  the  same  underestimate  of  stress  level  developed  at  long  times. 

(See  Figure  33). 
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COMPUTER    PROGRAM   "RELAX" 


FORTRAN  IV  G  LEVEl  1,  MOD  1 


0003 


/, 1040(2001 


0004 
0005 
0C06 
0007 
00C8 
0009 
0010 
0011 
0012 
0013 
001* 
0015 
0016 
0017 


0018 
0019 


0020 
0021 
0022 
0023 
002* 
0025 
0026 
0027 
0028 
0029 
0030 
0031 
0032 
0033 
003* 


0035 
0036 
0037 
0038 
0039 
00*0 
00*1 
00*2 
00*3 
00** 
00*5 
00*6 
00*7 
00*8 


0C*9 
0050 
0051 
0052 
0053 
005* 
0055 


0056 
0057 
0C58 
0C59 
00*0 


0061 
0062 
0063 
OCt* 
0065 
0066 
0067 
0068 
0069 
007C 
0071 
0072 
0073 


REAL  K,N,M,LABFL/*H 

DIMENSION  PLAST<5),STRI50I,T(200I,SIG(200I,TAI200),SIGA(200I,SIGB( 
12001, TB I  200  I, TH I  200). AREA (50 1 

READ  IN  THE  FOLLOWING  DATA 

tTITLE  »  TITLE  TO  BE  USED  ON  GRAPHICAL  OUTPUT 

PLAST  »  NAME  OF  MATERIAL  TESTED 

NDAT  -  NUMBER  OF  INITIAL  STRESS  CONDITIONS  FOR  WHICH 

THE  RELAXATION  CURVE  IS  TO  BE  INVESTIGATED. 
K.CB.O 
P,R,N,M  *  PARAMETERS  IN  THE  EMPIRICAL  CREEP  FQUATION 

FOR  THE  MATERIALS  TESTED 
STR  =  ARRAY  OF  INITIAL  STRESS  VALUES 
AREA  =  ARRAY  OF  CROSS-SECT inNAL  AREAS  FOR  THE  EXPERIMENTAL 
RUNS,  CORRESPONDING  TO  THE  INITIAL  STRESS  VALUES. 

READ  13,1  I  TITLE! 1 1 ,  1  -  1,121 
13  F0RMAM6A81 

READ  10, (PLASM  I  I,  I  -  1,*),  NOAT 

10  FORMATI*A*,I*l 

READ  11,  K,B,P,C,Q,N,R,M 

11  FORNAM3Fl5.5,5F6.*l 

READ  12.ISTRUI,  I  '    l.NDATI 

12  F0RMAM8F10.0) 
KCUE  =  50 
KCUE  =  0 

IF(KCUE.E0.50t   GO  TO  16 
[FUCUE.E0.01  GO  TO  1* 

1*  READ  12,(AREA(I),  I  ■  l.NDATI 
16  CONTINUE 

THE  00  LOOP  "00  102"  COMPUTES  THE  RELAXATION  CURVE  FROM 
THE  EMPIRICAL  CREEP  EQUATION  FOR  EACH  OF  THE  INITIAL 
STRESS  CONDITIONS. 

DO  102   L  =  l.NOAT 
S  *  STRIL) 

CALCULATE  THE  INITIAL  STRAIN, EO.   GENERATE  200  TIMES  WITH 
THF  LONGEST  ELAPSED  TIME  APPROXIMATELY  30  HOURS.   STORE 
THESE  TIMES  IN  THE  ARRAY  "T». 

EO  =  K*IS»«N1 
TINC  =  1./3600. 
T(l)  *  TINC 
DO  32   I  *  2,60 

32  MI)  =  Mt-1)  ♦  TINC 
TINC  =  10./360C. 

00  33   I  -  61,115 

33  T(  I  I  =  MI-ll  ♦  TINC 
TINC  -  1./30. 

DO  3*   I  *  116,1*1 
3*  Ml)  =  MI-ll  »  TINC 
TINC  «  0.5 
DO  35   I  =  1*2,199 

35  M  I  )  =  Tl I-ll  ♦  TINC 
S  =  0.95*S 

THE  DO  LOOP  "00  35"  USES  A  NEWTON-RAPHSON  ITERATION  TO 
SOLVE  THE  EMPIRICAL  CREEP  EQUATION  IMPLICITLY  FOR  STRESS 
AT  EACH  TIME  IN  THE  ARRAY  "T".   THE  INITIAL  STRAIN  IS  EO. 
THE  FIRST  ESTIMATE  OF  THE  RELAXED  STRESS  AT  ELAPSED  TIME 
Till  =  1  SECONO  IS  0.95*S.   THE  FIRST  ESTIMATES  AT  LATER 
TIMES  IS  THE  RELAXEO  STRESS  AT  THE  PREVIOUS  ELAPSED  TIME. 
STRESSES  COMPUTED  FOR  THE  TIMES  T(I)  ARE  STORED  IN  THF 
ARRAY  NAMED  "SIG". 

DO  36   I  =  1,199 

LC  =  0 

OIFF  =  *0000. 

PWR  i  C«(T( l)*»QI 

IFIPWR  -  50.1  82.82,81 

81  0  =  1.0 
GO  TO  83 

82  0  =  1.  -  EXP(-PWR) 

83  LC  =  LC  ♦  1 
SOLD  =  S 

FUNCT  =  -EO  ♦  K«(S*»NI  ♦  B« ( S»»R l*T (  I  I  »  P»(S**M|*n 

DERIV  =  N»K*(S**(N-1.) I»R»B*(S«*(R-1. ) l»T( I )*M*P»( S»*(M-1. I  1*0 

S  =  S  -  (FIINCT/DER  IVI 

OIFF  =  ABSISOLD  -SI 

IF  THE  STRESS  IS  NOT  FOUND  WITHIN  0.1  PSI  IN  TWENTY 
ITERATIONS,  ASSUME  CONVERGENCE  FAILURE. 

IF(LC.GT.20I  GO  TO  87 
IFIS.LE.O.OI  S  =  S/2. 
IFIDIFF  -  0.11  36,83,83 

87  PRINT  8B,  I.LC 

88  FORMAT! lOX,'  CONVERGENCE  FAILURE,  I 
S  =  0.75*STR(LI 

36  SIGI  II  =  S 


'.I2,'  LC  = 


,12, //I 


COMPUTE  THE  "FRACTION  OF  RELAXATION",  E3UAL  TO  RFLAXED 
STRESS  DIVIDED  BY  INITIAL  STRESS,  AND  STORE  IN  THF 
ARRAY  NAMED  "SIGA". 


DO  38   I  =  2,199 
TA(  I  I  =  MI-ll 
38  SIGA( II  =  SIGI I-ll /STRIL) 
TAUI  -    0.0 
SIGA(l)  =  1.0 


USE  SUBROUTINE  ORAW  TO  PLOT  THE  PREDICTED  RELAXATION  CURVE 
OF  "SIGA  VS.  ELAPSED  TIME". 


PRINT  *0,  LAST 
*0  FORMAT!//, 10X,'  LAST  =  ', 

DO  39   I  =  2,200 

TBI  I)  =  T(I-1)»3600. 

TH(  I)  =  MI-ll 
19  SIGB! I)  =  SIGI I-ll 

SIGBI 1)  =  STR1LI 

TBI  II  =  0.0 

THI11  =  0.0 

PRINT    ** 
**    F0RMAM1HII 

IFIKCUF.E0.50I       GO    TO    *8 
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C  COMPUTE  THE  PREDICTEO  RELAXED  10*0  BASED  ON  SPECIMEN  ARE* 

C  SO  TH»T  EXPERIMENTAL  LO*0-TlME  OATA  MAY  BF  COMPARED  WITH 

C  PREDICTED  LOAD-TIME  BEHAVIOR. 
C 

007*  DO  31   I  =  1.200 

0075  31  LOAD(l)  =  SI GB I  I  I *ARE A! L ) 

0076  PRINT  47,AREA(L) 

0077  47  FORMAT!///,  20X,  'THE  SPECIMEN  AREA  IS',F7.5,«  SOIIARF  INCHES',////! 
C 

C  PRINT  OUT  THE  PREDICTEO  STRESS  RELAXATION  CURVE  IN 

C  TABULAR  FORM. 

C 

0078  PRINT  45, (PLASTI II ,  I  =  1,4) 

0079  PRINT  49,  ( TH(  I ) , TBI  I ) , S IGBI I ) , LOAD  I  I > ,  I  =  1,2001 

0080  49  F0RMAT(1H0,20X,F10.6,13X,F8.0,13X,F8.1,13X,F8.1 I 

0081  GO  TO  103 

0082  48  CONTINUE 

0083  PRINT  45.  (PLASTII).  I  -  1,41 

0084  45  FORMAT! 10X, 'THE  CALCULATED  STRESS  RELAXATION  BEHAVIOR  FOR  ',4A4,' 
US  LISTED  BELOW:  ',//,20X,  'TIME,  HOURS'  ,  10X  ,  •  T  I  ME  ,  SFCONOS  •  ,  10X  ,  •  STS 
2ESS.PS! ' ,10X.'  LOAD',//) 


2ESS.PS! ',10X,'  LOAD',//) 
PRINT  46.ITHU  I  ,TB  I  I  I  ,  S  IGB  I  I  1,1  »  1.2001 
46  FORMAT!lH0,20X,F10.6,l3X,F8.0,13X,F8.1) 
103  CONTINUE 


0085  PRINT  46.ITHU  l,TBII  I  r  S  IGBI  I  1,1  »  1,200) 

0086 

0087 

0088  102  CONTINUE 


93 


TABULATED    OUTPUT    FROM   "RELAX' 


THE  SPECIMEN  AREA  ISO. 04992  SQUARE  INCHES 


THE  CALCULATED  STRESS  RELAXATION  BEHAVIOR  FOR 
TIME, HOURS  TIME, SECONDS 


NYLON       IS  LISTFn  BELOW: 
STRESS, PSI  LOAD 


0.0 

0.000278 

0.000556 

0.000833 

0.001111 

0.001389 

0.001667 

0.001944 

0.002222 

0.002500 

0.002778 

0.003056 

0.003333 

0.003611 

0.003889 

0.004167 

0.004444 

0.004722 

0.005000 

0.005278 

0.005556 

0.005833 

0.006111 

0.006389 

0.006667 

0.006944 

0.007222 

0.007500 

0.007778 

0.008056 

0.008333 

0.008611 

0.008889 

0.009167 

0.009444 

0.009722 

0.010000 

0.010278 

0.010556 

0.010833 

0. 011111 

0.011389 

0.011667 

0.011944 

0.012222 

0.012500 

0.012778 

0.013056 

0.013333 

0.013611 

0.013889 

0.014167 

0.014444 

0.014722 

0.015000 

0.015278 

0.015555 

0.015833 

0.016111 

0.016389 

0.016667 

0.019444 

0.022222 


0.0 

1. 

2. 

3. 

4. 

5. 

6. 

7. 

8. 

9. 
10. 
11. 
12. 
13. 
14. 
15. 
16. 
17. 
18. 
19. 
20. 
21. 
22. 
23. 
24. 
25. 
26. 
27. 
28. 
29. 
30. 
31. 

32. 
33. 

34. 

35. 
36. 
37. 
38. 
39. 
40. 
41. 
42. 
43. 
44. 
45. 
46. 
47. 
48. 
49. 
50. 
51. 
52. 
53. 
54. 
55. 
56. 
57. 
SB. 
59. 
60. 
70. 
80. 


2000.0 
1999.1 
1998.3 
1997.5 
1996.7 
1996.0 
1995.2 
1994.5 
1993.7 
1992.9 
1992.2 
1991. S 
1990.7 
1990.0 
1989.3 
1988.5 
1987.8 
1987.1 
1986.4 
1985.7 
1984.9 
1984.2 
1983.5 
1982.8 
1982.1 
1981.4 
1980.7 
1980.0 
1979.3 
1978.6 
1977.9 
1977.2 
1976.5 
1975.8 
1973.2 
1974.5 
1973.8 
1973.1 
1972.4 
1971.7 
1971.1 
1970.4 
1969.7 
1969.0 
1968.4 
1967.7 
1967.0 
1966.4 
1965.7 
1965.0 
1964.4 
1963.7 
1963.1 
1962.4 
1961.7 
1961.1 
1960.4 
1959.8 
1959.1 
1958.5 
1957.8 
1951.4 
1945.0 


99.8 
99.8 
99.7 
99.7 
99.7 
99.6 
99.6 
99.6 
99.5 
99.5 
99.4 
99.4 
99.4 
99.3 
99.3 
99.3 
99.2 
99.2 
99.2 
99.1 
99.1 
99.0 
99.0 
99.0 
98.9 
98.9 
98.9 
98.8 
98.8 
98.8 
98.7 
98.7 

98.7 
98.6 
98.6 
98.6 
98.5 
98.5 
98.5 
98.4 
98.4 
98.4 
98.3 
98.3 
98.3 
98.2 
98.2 
98.2 
98.1 
98.1 
98.1 
98.0 
98.0 
98.0 
97.9 
97.9 
97.9 
97.8 
97.8 
97.8 
97.7 
97.4 
97.1 
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0.025000  90.  1938.3  "(•.« 

0.027778  100.  1932.7  96.5 

0.030555  110.  192*-. 7  96.2 

0.033333  120.  1970.7  95. <) 

0.036111  130.  1914.8  "5.6 

0.038889  1*0.  1909.1  95.1 

0.0*1667  150.  1903.3  Q5.0 

0.0*****  160.  1897.7  94.7 

0.0*7222  170.  1892.2  9*.* 

0.050000  180.  1886.7  94. 7 

0.052778  190.  1881.2  93.9 

0.055555  200.  1875.9  93.6 

0.058333  210.  1870.6  93.4 

0.061111  220.  1865.3  93.1 

0.063889  230.  1860.2  92.9 

0.066667  2*0.  1855.1  92.6 

0.069***  .250.  1850.0  92.3 

0.072222  260.  18*5.0  92.1 

0.075000  270.  18*0.!  91.8 

0.077778  280.  1835.2  91.6 

0.080555  290.  1830.3  91.* 

0.083333  300.  1825.5  91.1 

0.086111  310.  1820.8  90.9 

0.088889  320.  1816.1  90.7 

0.091666  330.  18U.5  90.* 

0.09****  3*0.  1806. o  90.2 

0.097222  350.  1802.3  OO.O 

0.100000  360.  1797.8  89.7 

0.102777  370.  1793.*  89.5 

0.105555  380.  1789.0  89.3 

0.108333  390.  178*. 6  89.1 

0.111111  *00.  1730.2  88.9 

0.113838  *10.  1776.0  88.6 

0.116666  *20.  1771.7  8".* 

0.119**4  *30.  1767.5  88.2 

0.\222iZ  **0.  1763. 3  P8.0 

0.12*999  450.  1759.7  87.8 

0.127777  460.  1755.1  87.6 

0.130555  470.  1751.?  87.4 

0.133333  *80.  17*7.0  87.2 

0.136110  *90.  1743.0  87. n 

0.138888  500.  1739.1  86.8 

0.1*1666  510.  1735.'  86.6 

0.1*****  520.  1731.3  86.4 

0.1*7221  530.  1727.4  86.7 

0.149999  5*0.  1723.6  86.0 

0.152777  550.  1719.8  ss.8 

0.155555  560.  1716.1  8<;.7 

0.158332  570.  1712.3  as. 5 

0.161110  580.  1708.6  85.3 

0.163888  590.  1705.0  85.1 

0.166666  600.  1701.3  84.9 

0.1694*3  610.  1697.7  84.7 

0.202777  730.  1656.7  82.7 

0.236110  850.  1619.1  80.8 

C. 269443  970.  1584.6  79.1 

0.302777  1090.  1552.8  77.5 

0.336110  1210.  15?3.2  76.0 

0.369443  1330.  1495.7  74.7 

0.402777  1*50.  1*69.9  73.4 

0.*36110  1570.  1**5.8  77. 2 

0.*69*43  1690.  1427.1  71.  n 

0.502777  1810.  1*01.7  70.0 

0.536110  1930.  1381.5  69.0 

0.569*43  2050.  1362.*  68.0 

0.602776  2170.  13**. 3  67.1 

0.636110  2290.  1327.1  66.2 

0.669443  2410.  1310.7  65.4 

0.702776  2530.  1295.0  64.6 

0.736110  2650.  1280.1  67.9 

0.769443  2770.  1265.9  63.2 

0.802776  2890.  1252.2  62.5 

0.836110  3010.  1239.2  61.9 
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0.869443 

1130. 

0.902776 

3250. 

0.936109 

3370. 

0.969443 

3490. 

1.002775 

3610. 

1.036108 

3730. 

1.536108 

5530. 

2.036108 

7330. 

2.536108 

9130. 

3.036108 

10930. 

3.536108 

12730. 

4.036108 

14530. 

4.536108 

16330. 

5.036108 

18130. 

5.536108 

19930. 

6.03*108 

21730. 

4.536108 

23130. 

7.036108 

25330. 

7.536108 

27130. 

8.036108 

28930. 

8.336108 

30730. 

9.036108 

32530. 

9.536108 

34330. 

10.036108 

36130. 

10.536108 

37930. 

11.036108 

39730. 

11.536108 

41530. 

12.036108 

43330. 

12.536108 

45130. 

13.036108 

46930. 

13.536108 

48730. 

14.036108 

50530. 

14.536108 

52330. 

15.036108 

54130. 

15.536108 

55930. 

16.036102 

57730. 

16.536102 

59530. 

17.036102 

61330. 

17.536102 

63130. 

18.036102 

64930. 

18.536102 

66730. 

19.036102 

68530. 

19.536102 

70330. 

20.036102 

72130. 

20.536102 

73930. 

21.036102 

75730. 

21.536102 

77530. 

22.036102 

79330. 

22.536102 

81130. 

23.C36102 

82930. 

23.536102 

84730. 

24.036102 

86530. 

24.536102 

88330. 

25.036102 

90130. 

25.536102 

91930. 

26.036102 

93730. 

26.536102 

95530. 

27.036102 

97330. 

27.536102 

99130. 

28.036102 

100930. 

28.536102 

102730. 

29.036102 

104530. 

29.536102 

106330. 

30.036102 

108130. 

GkAPH  TITLED 

HfcRING.F.S., 
FRACTION  OF 

>  BOX  21,  STRESS 
INITIAL  STRESS 

RELAXATION  OF  NYLON 
VS.  TIME, HOURS 

1226.6 
1214.6 
1203.0 
1191.9 
1111.2 
1170.8 
1050.9 
972.3 
916.7 
875.3 
843.5 
818. 4 
798.2 
781.7 
768.2 
756.8 
747.3 
739.3 
732.4 
726.6 
721.5 
717.2 
713.5 
710.2 
707.4 
704.9 
702.7 
700.9 
699.2 
697.7 
696.4 
695.2 
694.2 
693.3 
692.5 
691.7 
691.1 
690.5 
689.9 
689.4 
689.0 
688.6 
688.2 
687.8 
687.5 
687.2 
686.9 
686.6 
686.  4 
686.1 
685.9 
685.7 
685.5 
685.2 
685.0 
684.8 
684.7 
684.5 
684.3 
684.1 
683.9 
683.8 
683.6 
683.4 


61.2 
60.6 
60.0 
•S9.5 
«9.0 
58.4 
52.5 
48.5 
45.8 
43.7 
42.1 
40.9 
39.8 
39.0 
38.3 
37.8 
37.3 
36.9 
36.6 
36.3 
36.0 
35.8 
35.6 
35.5 
35.3 
35.2 
35.1 
35.0 
34.9 
34.8 
34.8 
34.7 
34.7 
34.6 
34.6 
34.5 
34.5 
34.5 
34.4 
34.4 
34.4 
34.4 
34.4 
34.3 
34.3 
34.3 
34.3 
34.3 
34.3 
34.2 
34.2 
34.2 
34.2 
34.2 
34.2 
34.2 
34.2 
34.2 
34.2 
34.1 
34.1 
34.1 
34.1 
34.1 


HAS    BEEN    PLOTTED. 
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APPENDIX  E 
PREDICTION  OF  STRESS-STRAIN  CURVES  USING  THE  EMPIRICAL  CREEP  EQUATION 
A  simple  method  for  predicting  creep  curves  from  stress-strain  curves 
is  outlined  in  Figure  34.   As  indicated  in  Figure  34(a),  tensile  tests  are 
run  at  a  series  of  constant  strain  rates.   A  cross-plot,  Figure  34(b),  is 
constructed  for  reciprocal  strain  rate  vs.  strain  at  a  constant  stress. 
By  the  graphical  integration  shown  in  Figure  34(b),  times  corresponding  to 
certain  strains  may  be  found  for  a  given  stress  level,  and  creep  curves 
plotted  as  illustrated  in  Figure  34(c).   This  method  implies  the  existence 
of  an  equation  of  state  relationship  between  strain,  stress,  and  strain 
rate,  but  does  not  explicitly  state  such  an  equation.   A  similar  scheme 
was  proposed  by  Marin  and  Pao  [30] ,  who  obtained  good  agreement  from  test 
data  on  polymethylmethacrylate. 

In  this  study,  the  empirical  creep  equation  was  used  to  predict  stress- 
strain  curves  by  reversing  the  procedure  outlined  in  Figure  34.   From  the 
empirical  creep  equation, 


(E.l) 


e=  KaM  +  bcr*.  t    +■  Po-M(l-£"ct  ) 


The  slope  (strain  rate)  at  any  time  for  a  given  stress  may  be  found  by 


(E.2)    7T  =  uCT   +  YG    CQt   JL 


For  a  specified  strain  rate  and  a  given  stress,  equation  (E.2)  may  be 
solved  implicitly  for  time,  and  this  result  substituted  in  equation  (E.l) 
to  obtain  the  corresponding  strain.   A  Fortran  IV  computer  program,  called 
SIGEPS,  was  used  to  perform  these  calculations.   Output  from  SIGEPS  con- 
sisted of  stress-strain  data  for  a  material  for  a  given  strain  rate. 

97 


(a.) 


(b) 


% 


tiroa 


(c) 


Figure   34 
Prediction  of  Creep  Curves  from  Stress-Strain  Data 
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Program  SIGEPS  and  typical  output  is  given  on  the  following  pages.   Figure 
35  displays  this  output  along  with  experimental  tensile  test  data  obtained 
at  the  same  strain  rate.   The  same  characteristic  was  observed  for  most 
materials  tested;  the  predicted  tensile  strain  is  greater  than  the  observed 
strain  for  a  given  stress  level.   The  deviation  from  the  predicted  curve 
becomes  greater  at  higher  stress  levels.   In  the  tensile  test,  stress 
changes  with  time,  and  internal  structure  changes  probably  develop.   These 
changes  are  not  accounted  for  in  the  creep  equation  used  for  the  predic- 
tion.  In  addition,  it  should  be  noted  that  the  tensile  tests  took  place 
over  a  time  span  of  about  10  minutes  whereas  the  time  duration  of  the  creep 
test  was  100  hours.   Consequently  the  prediction  of  stress-strain  curves  is 
based  on  the  creep  curve  at  very  short  times.   As  indicated  in  the  tabulat- 
ed output  from  program  FITDATA  (See  Appendix  B) ,  it  is  at  these  short  times 
where  poorest  fit  of  equation  (E.l)  usually  occurred. 

Figure  36  illustrates  an  exception  to  the  previous  discussion.   The 
predicted  stress-strain  curves  for  polycarbonate  showed  excellent  agree- 
ment with  tensile  test  data.   It  is  possible  that  structure  factor  effects 
in  polycarbonate  are  negligible  at  the  stress  levels  encountered.   Further- 
more, output  from  FITDATA  for  polycarbonate  showed  good  agreement  between 
the  empirical  creep  equation  and  the  test  data  at  times  under  10  minutes. 

It  would  be  useful  to  further  examine  the  feasibility  of  this  method 
by  utilizing  an  empirical  creep  equation  specifically  designed  to  give 
good  data  fits  at  short  times. 
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20  ZO  40 

Strain  ,  in/in  x  lO3 


SO 


Figure  35 
Predicted  Stress-Strain  Curve  for  Nylon 
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Strain  ,  i*/'m  *  \03 


Figure   36 


Predicted  Stress-Strain  Curve  for  Polycarbonate 
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COMPUTER    PROGRAM    "SIGEPS" 


FORTRAN    IV    G    LEVEL     li    MOD    I 


MAIN 


DATE  =  68132 


8882 


MENS  ION  EDOT(50),E(5OI , T ( 501  . S IG( 50  I , FP ( 50 , 50  I , PL  AST ( 10 ) 


0005 
0006 
0007 
0008 
0009 
0010 
0011 
0012 


0013 

001* 
0015 
0016 


0017 
0018 
0019 
0020 
0021 


0022 
0023 
002* 
0025 
0026 


0027 
0028 

0029 
0030 
0031 
0032 
0033 
0034 
0035 
0036 
0037 
0038 
0039 
00*0 
00*1 

88*1 

00** 
00*5 
00*6 


enNi 


REAO    8,     (PLASTU  I,     I 
FORMAT?***) 


It*) 


READ  IN  THE  FOLLOWING  DATA: 

JPTS  *  NO.  OF  DIFFERENT  STRAIN  RATES  FOR  WHICH  STRESS- 
STRAIN  CURVES  ARE  TO  BE  COMPUTED 
NDAT  "  NUMBER  OF  DIFFERENT  STRESSES  AT  WHICH  THE  STRAIN 

IS  TO  BE  EVALUATED. 
A1,A2.A3.C,0  -  PARAMETERS  IN  THE  EMPIRICAL  CREEP 

EQUATION  FOR  A  GIVEN  STRESS 
EDOT  =  ARRAY  OF  STRAIN  RATES  FOR  WHICH  THE  STRESS  - 

STRAIN  CURVE  IS  TO  BE  DRAWN 
STR  =  STRESS  CORRESPONDING  TO  THE  CREFP  EQUATION 
PARAMETERS 

READ  10,  JPTS,  NOAT 

FORMATI2I5I 

READ  1*,(ED0T(II,  I  =  l.JPTSt 

READ  12,  A1,A2,A3,C,0 

FORMATI3E15.5,  2F10.01 

READ  13,  STR 

FORMAT(FIO.O) 

F0RMATI7F10.0I 

ESTIMATE    THE    TIME     AT    WHICH    A    STRAIN    RATF    FOOT! II     OCCURS    ON 
THE    CREEP    CURVE. 

GUESS    =    12. *A1)/IED0TI 11*60.1 

00    31       I    -    l.JPTS 

PRINT    86,     I 

FORMAT!//, 20X, •     I    ■    ',     12) 

USE    FUNCTION    SUBPROGRAM    "TIME"    TO    FIND    THE    TIME    AT    WHICH 
A    STRAIN    RATE    EDOTII)     OCCURS    FOR    THE    CREEP    CURVE    GIVEN    BY 
THE    EQUATION    WHEN    STRESS    -    STR.        THEN    USE    FUNCTION    SUB- 
PROGRAM   "STRAIN"    TO    COMPUTE    THE    STRAIN    FOR    THE    GIVEN 
STRESS. 

Till     =    TlME(GUESS.Al.A2,A3,Cf0,EDOT(I)l 
E(I)    =    STRAIN(A1,A2,A3,C,Q,T(I  II 
IFITI I I.EQ.10000.)       GO    TO    771 
EP(LC.I)    =    E<  I  ) 
SIG(LC)    «    STR 


CONTINUE    THF    PROCESS    FOR 
Tl  I) 


•NDAT"    DIFFERENT    STRESS    LEVELS. 


31  GUESS  = 
*0  LC  =  LC* 

IFILC  -  NOATI  9,9,23 
23  PRINT  2* 
2*  FORMAT(lHl) 

PRINT  OUT  THE  PREDICTED  STRESS-STRAIN  CURVE  IN  TABULAR 
FORM  FOR  EACH  STRAIN  RATE. 

PRINT  25,  (PLAST! II,  I  -  I,*) 

25  FORMATI10X,'  THE  STRESS-STRAIN  DATA  FOR • ,*A*, • FOR  VARIOUS  STRAIN  R 
1ATES  IS  LISTED  BELOW:',///) 

PRINT  26 

26  F0RMATI3X,'  STRESS',  26X,  '  STRAIN  RATES,  I N/ I N/MI N' , // I 
PRINT  27.(E00T(II,  I  =  l.JPTSI 

27  FORMAT! 12X.5F15.  6,//) 
DO  32   J  =     l.NOAT 

PRINT  28,  Sir.l  J|  ,<  EP(  J,  I  )  ,  I  x  l.JPTS) 

32  PRINT  29 

28  FORMAT!///, F10.0.5F15. 6) 

29  FORMAT! 1H0I 
GO  TO  *2 

771  II  =  I 

PRINT  772,  LCI 

772  FORMAT!///, 10X, '  CONVERGENCE  FAILURE  AT  LC  =  ',12,'  ,1  -     ',12,///! 
DO  773   I  »  II, JPTS 

773  EP(LC.I)  *    0.0 
GO  TO  *0 

*2  CONTINUE 
ENO 


02 


FORTRAN  IV  G  LEVEL  1.  MOD  1 


0002 
000? 
000* 
0005 
0006 
0007 
0008 
0009 


FUNCTION    STRAIN!  Al  ,A?,A3,  CO,  I) 


OATF    =    MH? 


FOR  A  GIVEN  T I HF  (I)  THIS  FUNCTION  SUBPROGRAM  USES  THF 
EMPIRICAL  CREEP  EQUATION  FOR  A  GIVFN  STRFSS  TO  COMPUTE 
THF    CORRESPONDING    VALUE    OF     STRAIN. 


POWER    =    C  *  I  7  **0» 
IFtPOWER    -     50.)    62,62, 61 

61  PROOT    =     1.0 
GO    TO    61 

62  PROOT    =     l.-EXP(-POWER) 

63  STRAIN    =    Al    ♦    A2»7 *A3»PR00T 
RETURN 

END 


FORTRAN 
0001 


IV  G  LEVEL  1,  MOO  1 


OATE  =  6R13? 


0002 
0003 
0004 
0005 
0006 
0007 
0008 
0009 
0010 
0011 
0012 
0013 
0014 
0015 
0016 

0017 


OOlfl 
0019 
0020 
0021 
0022 
0023 
0024 


FUNCTION  TIME! X, Al , A2 , A3 ,C ,0,E0) 


THIS  FUNCTION  SUBPROGRAM  COMPUTES  THF  T I  MP  (X)  AT  WHICH  A 
GIVEN  STRAIN  RATE  (ED)  OCCURS,  BASED  ON  THE  CREEP  PARAM- 
ETERS A1,A2,A3,C,Q  FOR  A  GIVEN  STRESS.   THE  FUNCTION 
CALLED  "FUNCT"  IS  THE  DERIVATIVE  OF  THE  CREEP  EQUATION, 
OE/DT,  MINUS  THE  GIVEN  STRAIN  RATE, ED.   THF  SUBPROGRAM 
USES  A  NEWTON-RAPHSON  ITERATION  TO  FIND  THE  UNKNOWN  TIMF 
(X),  WHICH  IS  THE  ROOT  OF  THE  EOUAT I ON-FUNCT  =  0". 
"DERIV"  IS  USED  IN  THE  N/R  ITERATION  AND  IS  THE  DERIVATIVE 
OF  "FUNCT"  TAKEN  ANALYTICALLY. 

L  =  0 

R  -  Q  -  1. 

90  XOLD  -  X 

PRINT  07, L, XOLD 
87  FORMAT! 25X.'  L  =  '  ,12, 

PWR  =  C»IX**0) 

1FIPWR    -    50.)     80,80,82 
80    PROD    -    A3»EXPI-PWRI 

FUNCT    =    A2    -    ED    *    (C*0»PROD* I X»*R I  I 


XOLD    =    ■ ,F15.9I 


DERIV    =    PROD*C*0»(R»( X»*(R-1. ) ) 
X    =    X    -     IFUNCT/DER IVI 
IFIX.LT.O.OI        X    =X0LD/2. 
DIFF    =    ABSIX    -    XOLD) 
PRINT    93, OIFF, FUNCT, OERIV, PWR 
93    F0RMATI25X,1    DIFF    =     ',F15.9,' 
1,'  PWR    =     '.F10.6) 

L    =    LM 


C*0*(X»*(?.*R) ) ) 


,F!5.9, •     DERIV    = 


RET 
END 


L.GT. 
DIFF 
E  =  ) 
URN 
E  =  ] 
URN 


IF    more    THAN    20    ITERATIONS    ARE    USED    TO    FIND    THE 
TIME     (X),     ASSUME    CONVERGENCE    FAILURE. 

201       GO    TO    82 

-     .000011     81,81,90 


TABULATED     OUTPUT    FROM    "SIGEPS" 

THE     STRESS-STRAIN    DATA    FOR       POLYCARBONATE    FOR    VARIOUS    STRAIN    RATES     IS    LISTFD    BELOW: 

STRESS  STRAIN    RATES,     IN/IN/MIN 

0.013040 

2000.  0.0060»7 

2500.  0.007790 
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APPENDIX  F 
THE  BOX  DISTRIBUTION  FUNCTION  FOR  RETARDATION  SPECTRA 

The  use  of  a  continuous  distribution  of  retardation  times  to  character- 
ize delayed  elastic  response  was  discussed  in  Section  6.   Not  only  do 
theoretical  considerations  favor  this  approach,  but  mathematical  solutions 
to  viscoelastic  problems  may  be  considerably  simplified,  particularly  if 
the  functions  representing  the  spectra  can  be  handled  by  transform  techni- 
ques.  As  an  example,  an  analytic  expression  for  creep  behavior  will  be 
derived  here,  using  a  box  distribution  of  retardation  times  as  illustrated 
in  Figure  37 . 

An  expression  for  creep  strain  corresponding  to  the  four-element  model 
of  Figure  10,  page  25  ,  may  be  written: 


(F.D       e  -    A,  +    Alt     +  A?(\-JLbt) 


where  the  constant  D  is  the  rate  constant  equal  to  the  reciprocal  of  the 
characteristic  relaxation  time  T    ,    and  coefficients  A..  ,  A„,  and  A„,  de- 
pend on  the  material  and  the  applied  stress.   For  a  material  with  many  dis- 
creet retardation  times,  the  expression  for  creep  strain  given  by  equation 
(F.l)  becomes: 


-bit 


(f.2)        g  =    A,   +  Azt    +  A*Tf  a.i(\-JL~  l  ) 


where  the  constant   t>t   is  the  rate  constant  associated  with  the   L 


+h 


retardation  time   "77*  ,  and  the  factors  (X,i     are  weighting  constants  in- 
dicating  the  relative  contribution  of  the    L —   mechanism  to  the  delayed 
strain.   Since  the  factors  GLi    are  chosen  so  that  1E.Q-L  ='  1  ,  equation 
(F.2)  may  be  rearranged  as: 
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Figure  37 
Box  Distribution  of  Rate  Constants 


<j(b) 


Figure  38 
Wedge  Distribution  of  Rate  Constants 
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€4.— 


(f.3)     e  =   A,    +  Att    +■    A3(\-  Ucrb;t) 

If  the  delayed  elastic  response  (  £.d   )  has  associated  with  it  a  continu- 
ous spectrum  of  retardation  times,  then  the  factors  CL\    may  be  represent- 
ed by  a  rate  constant  distribution  function,  ^L(  b)   ,  such  that: 

(F.4)     J   <j(b)  cCb   -  \ 
For  this  continuous  distribution  of  retardation  times  the  creep  equation 


(F.3)   becomes: 


(F.5) 


bt 


£  =    A,  +  Att  +  A*(l- j^b)^  tftb) 


Equation  (F.5)  now  contains  the  Laplace  transformation  of  the  rate  con- 
stant distribution  function   9,(t>)  »  anc*  may  be  rewritten  as: 

or.6)       e  =   A,  +  A,t  +  A,(|-  £^(b)]) 

where   Xv^C^M   represents  the  Laplace  transform  of  the  function. 

The  transform   i-\^(bM  can  be  evaluated  readily  in  the  case  of 
a  box  distribution  of  rate  constants.   Referring  to  the  notation  used 
for  the  box  distribution  illustrated  in  Figure  37, 


ft*"} »  L^£btdh 


In  equation  (F. 7) , constants  C~  and  C-  are  the  maximum  and  minimum  rate 
constants,  respectively,  of  the  box  distribution.   These  constants  may 

106 


be  evaluated  for  a  given  set  of  creep  data  as  shown  in  Figure  19. 

In  equations  (F.l)  and  (F.6),  the  coefficient  A„  is  intended  to 
represent  the  delayed  strain  after  elastic  equilibrium  is  achieved,  i.e., 
after  infinite  time.   The  delayed  strain  at  time  zero  is  zero.   Therefore, 
the  expression  for   dv_\^  "^ '  V    given  by  equation  (F.7)  must  be  normalized 
so  that: 

(F.8)    t{<j(bV}  =  O   ,   t=0 


■—   oo 


This  may  be  accomplished  by  expanding  the  exponential  terms  of  equation 
(F.7)  in  a  MacLaurin  series: 

(F'10)    fttftfW}-   ct-c, 

Equation  (F.7)  may  now  be  normalized  to  get: 

(F.u)  lw  =  fciW-t U"c,t - *"] 

By  substitution  of  the  time  limits  L  -  O   and  "C  —  °°   into  equation 
(F.ll)  it  is  apparent  that  the  conditions  of  equations  (F.8)  and  (F.9) 
are  satisfied.   Substitution  of  equation  (F.ll)  into  equation  (F.6)  gives 
the  analytic  expression  for  creep  strain  when  the  retardation  spectrum  is 
a  box  distribution  of  rate  constants. 


(F.i2)    e  =  A,  +  A2t    +  A 


(cx-c,yt 
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Similar  derivations  may  be  performed  whenever  the  retardation  spectrum 
can  be  represented  by  a  sum  of  functions  whose  transforms  can  be  found. 
Unfortunately,  the  resulting  expression  is  quite  complex  for  all  but  the 
simplest  distribution  functions.   For  example,  the  wedge  distribution  func- 
tion, Figure  38,  similar  to  the  theoretical  model  of  Alfrey  discussed  in 
Section  6,  leads  to  the  expression: 

(,.13)    £=A,  +  A2t  +  A3[l-  ^t\K-K.^r\ 

In  the  equations  for  creep  strain  developed  in  this  Appendix,  the 
delayed  elastic  response,   €iol  ,  represents  a  portion  of  the  total  strain. 
In  equation  (F.l),  for  example, 

(F.14)   Ga  «  A^O  --£    ) 
The  value  of  the  strain  €.<£  varies  from  zero  at  "t  =■  O   to  a  value  A 
at  "t  -  °°  -      It  is  often  convenient  to  normalize  the  delayed  strain.   The 
normalized  form,  called  the  delayed  strain  fraction,  is  defined  by: 

(f.15)    £  =  eyA^ 

The  delayed  strain  fraction  varies  from  £  =  0  at  ~t  -  O   to  £  —  I   at 
L    =  °°   ,  and  its  value  at  any  time  L   >  O     represents  the  progress 
toward  elastic  equilibrium.   For  equation  (F.l),  the  delayed  strain  frac- 
tion is: 

(F.16)    £.  -  \  -  SL 

A  better  insight  into  the  nature  of  the  delayed  elastic  response  is 
achieved  by  plotting  delayed  strain  fraction  versus  the  logarithm  of  time. 
Equation  (F.16)  is  plotted  in  this  manner  in  Figure  39.   The  inflection 
point  of  this  sigmoidal  shaped  curve  is  reached  at  a  time  t  =  /[^  ,  the 
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characteristic  retardation  time.   The  delayed  strain  fraction  at  "t  =  /b 
is  £  =  \~~   yj?_   ~   0.632.   The  slope  of  this  curve  changes  very  gradually 
in  the  range   t  =  0.3  to  £=0.8.   In  this  region,  the  shape  of  the  curve 
can  be  approximated  by  a  straight  line  drawn  with  a  slope  corresponding  to 
<A£./cL(t»i.b)         evaluated  at   £  =  0.5.   This  is  shown  by  line  segment  A-A 
in  Figure  39. 

For  a  box  distribution  of  rate  constants,  the  delayed  strain  fraction 
versus  log  time  curve  (Figure  40)  is  also  sigmoidal  in  shape,  with  even 
less  drastic  slope  changes  in  the  region  around   <£  =0.5.   In  order  to 
accurately  depict  the  delayed  elastic  response,  the  sigmoidal  curve  must 
fit  the  data  closely  at  a  time  t4  corresponding  to   £  =  0.5.   Otherwise, 
large  errors  in  predicted  strain  occur  over  most  of  the  time  range.   Fig- 
ure 41(a)  illustrates  a  sigmoidal  curve  exhibiting  a  slope  mismatch  with  a 
data  plot.   In  Figure  41(b),  the  slopes  of  the  data  plot  and  the  sigmoidal 
function  are  the  same,  but  the  responses  do  not  occur  at  the  same  time. 
It  is  apparent  from  these  figures  that  the  error  in  predicted  delayed 
strain  fraction  will  be  large  if  there  is  a  slope  and/or  time  mismatch  at 
approximately   £.  =0.5. 

Plots  of  £.  versus  /C4^L   were  made  for  the  creep  data  obtained  in 
this  investigation.   For  different  viscoelastic  materials  tested,  these 
plots  exhibited  remarkable  similarities.   The  delayed  strain  fraction 
reached  a  value   £L  =  0.5  at  a  time   t  *  1  hour,  and  the  slope  of  the 
curves  at  £.  =  0.5  was  d£/cL(hx)   «*  0.15.   The  same  values  were  observed 
for  PP,  PC,  NY,  and  PE  and  were  independent  of  stress.   Therefore,  the 
same  function  may  be  used  to  characterize  the  delayed  elastic  response 
of  each  of  these  materials.   A  question  arises  as  to  whether  or  not  a  func- 
tion such  as  the  box  distribution  leads  to  a  suitable  representation.   To 
text  the  applicability  of  the  box  distribution,  it  is  necessary  to  investi- 
gate the  characteristics  of  the  resulting  sigmoidal  curve  for  the  limiting 
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Figure  39 
Delayed  Strain  Fraction  vs.  Ln  Time  for  a  Single  Retardation  Time 


Figure  40 
Sigmoidal  Curve  for  a  Box  Distribution 
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Figure  41 
Unsuitable  Representations  of  Delayed  Elastic  Response 
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cases  of  the  box  distribution  shown  in  Figure  42.   This  analysis  is  given 
below: 

If  the  limits  of  the  box  distribution  are  extended  to  cover  many  de- 
cades of  time,  as  shown  schematically  in  Figure  42(a), 

£ Clt  «  XC,t       ^r   t  >  O 
The  expression  for  the  delayed  strain  fraction  (from  equation  (F.12))  is 
given  by: 

(F.17)    C  -   \     \_      (c.-C.Vt 

-C   t 
The  terms  J£,      *•       and  C.  may  be  neglected  so  that, 


(F.i8)     a.  =    i  - 


■c,t 


Ctt 


The  constant  C1  is  very  small  so  that, 

The  delayed  strain  fraction  for  the  extended  box  distribution  may  then  be 
approximated  by: 

(F.19)    £  =   I  -   Clt     i   Cl>>C' 
This  curve  is  illustrated  in  Figure  43(a)..   The  approximation  is  obvious- 
ly in  error  for  times   L  ^  /Cf_   >  but  the  constant  C„  is  large  and  the 
equation  (F.19)  is  valid  in  the  time  region  of  interest. 

For  the  case  shown  in  Figure  42(b)  where  the  limits  converge  and  the 
box  becomes  narrow,  the  equation  for  delayed  strain  may  be  approximated 
as  follows:   Denoting  the  bracketed  term  of  equation  (F.17)  by  ZL   ,  and 
expanding  the  exponentials, 

,    (&tV    (gip   (Clfc?  fctV   CczOVM)* 
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(P.20)    I  ^  <^=£iL   -  CCiSrt  ..t  ^  £*-c?)    t     ... 

•16. -co    2.CC-C0     3TTrx^c7) 

<i~     torn  of  this  expansion  in  represented  by: 
C,A1  -C.^   ,  M-J 


The  va1    >     Ls  term  for  the  case  of  conver  ing  limits  (.  C?   -      --   C J 
may  be  determine  J  by  considering  C9  and  C,  to  differ  by  a  small  quantity  f 


then 


c2  —  >  cl    ,   c2  =  c. 


a,   =   (c1  +  e  ) 


Using  a  binomial  expansion  and  truncating  after  the  second  term, 
therefore, 


Substituting  equation  (F.21)  for  the  /*.  ten-,)  of  equation  (F.20)  and 
evaluating  the  terms  for  /YL~   1,2,3 .the  expansion  becomes: 

j.       •  ,  z.  ^-' 

C.t"     c^t'         C     .*-» 

Enuation  (F.2?)  represents  the  expansion  of  j£        in  a  MacLaurin 

series,  so  the  desired  approximation  has  been  found: 

(F.23)       £ ■  *  '  ~  £ 

-  c»t 

or  f  s  |  -  £  for  C2  • — *  C 
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Equation  (F.23)  would  result  from  a  "single  spike"  distribution,  corres- 
ponding to  a  single  characteristic  retardation  time.   The  analysis  veri- 
fies the  expected  result  that  equation  (F.17)  should  reduce  to  equation 
(F.23)  as  the  width  of  the  box  distribution  becomes  vanishingly  small. 
This  approximation  is  sketched  in  Figure  43(b). 

In  order  to  be  used  for  an  accurate  representation  of  the  data  in  this 
investigation,  the  slope  of  the  sigmoidal  plot  of  the  box  distribution  must 
take  on  the  value : 


JLZ_ 

(F-24)  tKUt) 


rz        0.\S 


By   the   approximation   for  the   limiting  case   of  C?   ^>    C. ,    given  by  equa- 
tion  (F.19), 

1 

£=  \  -  5"t 

for  d  «0.S  ,    t=  I 
Ct  =   Z 

Jit    "    C,"tx     )      cLfaJ:)  ~  di 
oU  JL 


<F-25>    ie_ 

cL(fa±) 


£  =  0.5         —        0»S 

t-\ 


For  the  case  of  narrowing  box  width,  (C_  — *>  C..),  given  by  equation 
(F.23):  t 


At 


^T  =  Ci£ 


-<*t 


for  £=0.5,t-l 

£•*.*=  0.7 
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ct  >»  c. 


fft) 
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b  — 


b  — 


(a) 
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Figure  42 
Limiting  Cases  of  the  Box  Distribution 
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Figure  43 

Approximations  of  the  Response  Function  for  the 
Limiting  Cases  of  the  Box  Distribution 
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gig-    -    r    tJL 


-ctt 


JL(JU±) 


£=o.S   «:  0.35 


It  may  be  demonstrated  by  numerical  methods  that  equation  (F.26)  represents 
the  minimum  value  of  the  slope  d.S./eL[X^t)  for  any  choice  of  the 

constants  C.  and  C-,  provided  C~   >   C  .   The  requirements  of  equation 
(F.24)  can  not  be  met,  since  the  observed  slope  of  0.15  does  not  fall  in 
the  range  of  0.35  to  0.5.   Therefore  the  box  distribution  is  an  unsatis- 
factory choice  for  the  retardation  spectrum  of  the  materials  tested  in 
this  experiment. 

A  similar  analysis  of  the  wedge  distribution  (Figure  38)  shows  a 
slope  range  0.35  ^lci£/cL(&d-)  k=o.s)  ^  0.4  for  all  choices  of  K.  and  K_, 
indicating  that  this  rate  constant  distribution  is  also  unsatisfactory. 
Since  the  limiting  case  of  the  box  distribution  given  by  equation  (F.23) 
corresponds  to  the  response  of  a  single  Voigt  Element  (Figure  5),  the  four 
parameter  creep  model  (Equation  (F.l))  gives  a  poor  representation  of  de- 
layed elastic  strain. 

As  discussed  in  Section  6,  good  correlations  can  be  obtained  using 
a  model  with  several  discreet  retardation  times.   This  "multi-spike  distri- 
bution" produces  a  plot  of  delayed  strain  fraction  versus  log  time  with 
inflections  at  each  retardation  time.   These  retardation  times  may  be  chosen 
so  that  the  response  function  approximates  the  creep  data  quite  closely. 
Figure  44  illustrates  the  response  function  resulting  from  a  three-spike 
distribution  superimposed  on  a  typical  plot  of  creep  data.   The  response 
function  from  a  box  distribution  is  also  shown.   The  lower  limit,  C  ,  of 
the  box  is  equal  to  the  slow  rate  constant  Dj  of  the  three-spike  distri- 
bution, and  the  upper  limit  C-  equals  the  fast  rate  constant  D3  . 
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These  values  for  the  box  distribution  function  were  chosen  to  illustrate 
the  dominant  role  of  the  fast  rate  constants  in  governing  the  delayed 
elastic  response.   The  constants  C.  and  C     could  be  changed  to  shift  the 
sigmoidal  curve  toward  longer  times,  but  the  slope  cannot  be  adjusted  to 
produce  a  good  fit  with  the  data. 

Monomodal  distribution  functions  characterized  by  a  single  maximum 
can  always  be  approximated  by  a  sum  of  box-type  distributions  as  shown  in 
Figure  45.   The  limits  of  each  succeeding  box  fall  within  the  range  of  the 
first  box.   Since  the  slope  of  the  delayed  strain  fraction  versus  log 
time  plot  resulting  from  such  a  distribution  is  restricted  in  accordance 
with  equations  (F.25)  and  (F.26),  it  is  apparent  that  the  retardation 
spectrum  of  the  plastics  tested  is  not  monomodal. 

Bimodal  (or  polymodal)  spectra  produce  a  series  of  inflections  in  the 
plot  of   £_  vs.  -^-vtt  ,  which  results  in  a  lower  overall  slope,  i.e.,  the 
delayed  elastic  response  occurs  over  many  decades  of  time.  The  three- 
spike  distribution  of  Figure  44  is  an  example  of  a  polymodal  retardation 
spectrum. 

The  delayed  elastic  response  of  the  materials  tested  covered  about 
five  decades  of  time.   This  response  can  be  characterized  by  a  polymodal 
retardation  spectrum  such  as  illustrated  in  Figure  46,  where  the  peaks 
imply  the  presence  of  several  different  molecular  mechanisms  whose  time 
constants  differ  by  orders  of  magnitude.   This  model  is  in  accord  with 
Alfrey's  concept  (See  Section  6)  of  local  bond  straightening  and  long- 
range  uncoiling  of  the  polymer  chains. 
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Figure  45 
Monomodal  Retardation  Spectrum  Approximated  By  Box  Distributions 


Figure  46 
A  Polymodal  Distribution  Function 
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APPENDIX  G 
EQUIPMENT 

The  creep  machines  used  in  this  investigation  apply  constant  dead- 
weight loading  through  calibrated  lever  arms.   (Figure  47).   Eleven  identi- 
cal units  are  installed  in  one  test  stand.   Photographs  of  the  laboratory 
and  apparatus  are  shown  in  Figures  48  through  50.   Specimens  were  gripped 
by  steel  pins  placed  through  the  5/8  inch  diameter  hole  in  each  end,  and 
uniaxial  loading  was  maintained  through  a  series  of  universal  joints 
above  and  below  the  grips.   Turnbuckles  and  an  adjustable  counterbalance 
weight  permitted  final  balancing  of  the  lever  arm  prior  to  load  applica- 
tion. 

Studies  of  viscoelastic  behavior  of  materials  require  large  amounts 
of  data  obtained  over  long  time  periods  of  testing.   Consequently,  facili- 
ties for  automatic  recording  of  creep  and  relaxation  data  are  highly  de- 
sirable.  Design  of  such  instrumentation  was  started  in  conjunction  with 
this  project  and  is  presently  being  installed.   Special  clip-on  strain 
gage  fixtures  were  fabricated  to  permit  measurement  of  the  large  strains 
encountered  in  mechanical  testing  of  viscoelastic  materials.   These  clip 
gauges,  shown  in  Figures  51  and  52  consist  of  phosphor  bronze  strips  with 
BLH  SR-4  Type  A-5-1  strain  gauges  mounted  on  each  side.   Tests  on  phosphor 
bronze  indicate  that  this  material  has  remarkable  linearity  over  the  en- 
tire elastic  range  with  no  observable  hysteresis  effects.   The  strips  are 
connected  to  aluminum  legs  which  make  a  knife-edge  contact  with  the  speci- 
men.  This  design  has  two  useful  features: 

(1)  The  legs  are  detachable  so  that  different  leg  lengths  may 
be  used  depending  on  the  strain  sensitivity  desired. 

(2)  The  mode  of  attachment  to  the  specimen  is  such  that  this 

electromechanical  device  may  be  used  for  specimens  of  various 

widths  and  thicknesses. 
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The  four  strain  gauges  are  connected  in  a  Wheatstone  bridge  arrange- 
ment shown  in  Figure  53.   A  potentiometer  in  the  circuit  permits  a  null 
adjustment  to  be  made.   This  bridge  circuitry  has  several  advantages:  [26] 

(1)  Automatic  moment  compensation  (in  the  event  there  are  devia- 
tions from  axial  loading. 

(2)  Signal  strength  increased  by  a  factor  of  four,  (over  that 
of  a  single  gauge  mounted  on  the  specimen) 

Terminal  boards  have  been  set  up  on  the  eleven-unit  creep  machine  to 
receive  the  leads  from  the  extensometer  circuit  for  each  test  specimen. 
Bridge  circuits,  each  with  its  own  balancing  potentiometer,  are  placed 
in  parallel  across  a  D.C.  power  supply,  with  leads  from  each  bridge  enter- 
ing master  control  panels.  From  these  panels,  the  output  signal  from 
each  extensometer  may  be  fed  into  a  multi-channel  recorder,  permitting 
testing  and  graphical  data  output  for  up  to  eleven  specimens  simultaneous- 
ly.  A  stress  analysis  of  the  clip  gauges,  using  linear  elastic  theory, 
is  given  in  Figures  54  and  55. 

It  is  desirable  to  have  the  output  of  electromechanical  extensometers 
converted  directly  into  digitalized  form  because  of  the  usefulness  of  the 
digital  computer  in  data  analysis.   Preliminary  design  of  such  instru- 
mentation has  begun  at  NPGS,  Monterey. 
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Figure  47 
Lever  Arm  Loading  Assembly 
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Figure  48 

Specimen  Grips  and  Clamps 
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Figure  49 

Eleven  Unit  Creep  Testing  Apparatus 


Figure  50 
Tripolitis  Extensometer  Attached  to  Specimen 
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Figure  53 
Bridge  Circuitry  for  Extensometer 
125 


1 


!,      1 

5R-4  gaga*  -/       /               ; 

' 

' 

/                                           ^ 

h     ; 

Z_  phosphor          x 

' 

bror»z.a           ^ 

J                                   5tr\p               J 

/* 

s 

* 

s 

s 

s 

\ 

V 

aluminum 


Deformed 
Gaqa 


V>b 


f 


Strip  Cross 
±       faction 


* 


"T 


Figure  54 
Dimensioned  Sketch  for  Clip  Gauge  Stress  Analysis 
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CLIP  GAUGE  STRESS  ANALYSIS 


•  Assume  no  bending  in  aluminum  legs 

•  To  find  €.  in  gauges  for  a  given  specimen  extension,  S 

pe  =  ~    ,  e-Ai^0\Th) 

from  the  flexure  formula,    £  —  ^/P 

for  small  angles    0  ,       AJuZ   (^jrj   «   ~i 

for    cU.Olo",    b«Z",    L-l" 
£  =  .00115  S 

for  a   specimen  extension  of    $  =  |" 
£  =  12.50/tincbw/inch 

or  about  one -ha If  the  maximum  strain 
capability  of  the  BLH  type  a- 5-1  gages 

maximum  stress  in  the  bronze  strip: 
/r-   Zc  €*.00»-LS 

u  ~  tfe  ;      E=  J6*\0*pa\ 

(S  -  ZOpOO  p*\ 

or  only  a  small  fraction  of  the  yield 

strength  of  phosphor  bronze  (  Oy  =  112,000  psi) 

at  an  extension  of  S—  I  "  ,  error  in  the  strain 
due  to  the  small  angle  assumption  is: 

VtW=  'A*  =•«  ,  AM? (&l  =.M74 

,  Error  =   A»-1*T4  «  p/ 
.15        A  '© 


Figure  55 
Theoretical  Stress  Analysis  of  Clip  Gauges 
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